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/  ABSTRACT 


'y  The  external  nucleate  pool-boiling  heat-transfer  coeffi¬ 
cient  of  a  horizontal  smooth  copper  tube  in  R-1 14-oil 
mixtures  (0  to  10  percent  oil)  was  measured  for  heat  fluxes 
from  1  to  100  kW/m5"  at  two  different  saturation  temperatures 
(-2.2  «C  and  6.7  °C)  .  A  copper-nickel  tube  coated  with  the 
Onion  Carbide  "High  Flux"  coating  was  similarly  tested.  The 
High  Flux  coating  was  found  to  improve  the  heat-transfer 
coefficient  by  at  least  a  factor  of  7  in  oil- free  R-114.  Oil 
resulted  in  about  a  20  percent  reduction  of  the  heat- 
transfer  coefficient  of  the  High  Flux  surface  at  heat  fluxes 
less  than  30  kW/m2  and  up  to  an  80  percent  reduction  at  heat 
fluxes  above  30  kW/m2  with  greater  than  6  percent  oil. 
Onder  all  conditions,  the  High  Flux  coated  tube  outperformed 
the  smooth  copper  tube.  (  r  ■  ■  -  .  u  ; 
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I.  introduction 


1.  BACKGROUND 

The  U.S.  Navy  currently  uses  refrigerant  R-114  in 
centrifugal  chilled-water  air-conditioning  plants  aboard 
submarines  and  surface  ships.  The  Navy  hopes  to  reduce  the 
sire  of  these  units  and  increase  their  performance  by  using 
enhanced  evaporator  and  condenser  surfaces.  An  experiment 
by  Aral  et  al.  [Ref.  1]  produced  a  prototype  200-Ton  R-12 
centrifugal  water  chiller  that  was  28  percent  shorter  in 
length  and  had  a  50-70  percent  improvement  in  the  overall 
heat-transfer  coefficient  by  employing  the  enhanced  surface 
"Thermoexcel  E"  made  by  the  Hitachi  Company.  Comparisons  of 
various  enhanced  commercial  tubes  by  lilmaz  and  West  water 
[Ref.  2],  Marto  and  Lepere  [Ref.  3],  and  Carnavos  [Ref.  4] 
for  various  refrigerants  other  than  R-114  indicated  that  a 
porous-coating-enhanced  surface,  such  as  Union  Carbide's 
"High  Flux,"  will  exhibit  the  best  boiling  heat- transfer 
performance  in  a  pure  refrigerant. 

The  High  Flux  surface  (see  Figure  1.1)  consists  of  a 
sintered  metallic  matrix  bonded  to  a  metallic  substrate. 
The  surface  is  produced  by  coating  a  smooth  tube  with  a 
binder- solvent  mixture  and  then  applying  a  mixture  of  metal 
and  braze  alloy  powder;  the  tube  is  placed  in  a  furnace  to 
evaporate  the  solvent,  binder,  and  melt  braze  alloy  thus 
forming  a  porous  structure  having  multiple  reentrant  cavi¬ 
ties  to  enhance  nucleation.  [Ref.  5] 

Since  the  High  Flux  surface  will  be  employed  in  a 
refrigeration  unit  using  an  oil-lubricated,  hermetically- 
sealed,  compressor,  some  amount  of  oil  is  always  present  in 
the  evaporator.  Studies  by  Henrici  and  Hesse  [Ref.  6]  for 


VAPOR  BUBBLE  TRAPPED 
IN  ACTIVE  SITE 


smooth  tubes  and  Stephan  £  Ref-  7]  for  the  Gewa-T  surface 
(manufactured  by  the  Wieland  Company)  in  R-114-oil  mixtures 
indicate  that  the  heat-transfer  coefficient  of  enhanced 
surfaces  can  be  significantly  altered  by  oil. 

Experimental  data  showing  the  effect  of  oil  concentra¬ 
tion  on  the  heat-transfer  coefficient  of  the  High  Flux 
surface  in  R-114  are  lacking,  thus  motivating  the  present 
investigation.  This  investigation  was  funded  by  the  David 
W.  Taylor  Naval  Ship  Research  and  Development  Center. 
Details  of  the  experimental  apparatus  are  described  by 
Karasabun  [Ref.  8].  The  smooth  copper  tubes  were  supplied 
by  the  Nieland  Company.  The  High  Flux  coated  copper-nickel 
tubes  were  supplied  by  the  Onion  Carbide  Corporation. 

B.  TBES1S  OBJECTIVES 

The  objectives  of  this  thesis  are: 

1.  Take  boiling  data  on  a  smooth  tube  in  R-114  with  and 
without  oil  for  comparison  with  the  data  of  other 
researchers,  and  to  provide  baseline  data  for  evalu¬ 
ating  the  boiling  performance  of  the  High  Flux  tube. 

2.  Take  boiling  data  on  a  High  Flux  tube  for  various  oil 
concentrations  (0  to  10  percent  by  mass). 

3.  Study  the  effect  of  saturation  temperature  on  the 
R-114  boiling  behavior. 

4.  Attempt  to  sample  oil  locally  in  the  near  vicinity  of 
a  tube  to  investigate  the  possibility  of  an  oil 
concentration  gradient  around  the  tube  during 
operation. 


II-  HEYIEi  OF  BEF8IGE BAHT-OIL  Mil TO HE  BEBAVIOB 

A.  BUCL EATE  BOILIHG  07  BEFBIGEBAHT-OIL  MIXTURES  FBOM  SMOOTH 

TUBES 

In  1963,  Stephan  [Bef.  9]  published  a  milestone  paper  on 
the  influence  of  oil  on  the  boiling  heat  transfer  of  B-12. 
The  effects  he  noted  have  been  observed  in  most  refrige¬ 
rants,  including  R-114.  In  1972,  Henrici  and  Hesse  [Bef.  6] 
updated  Stephan* s  work  for  R-1 14-oil  mixtures  boiling  from  a 
smooth  copper  tube.  Figures  2.1  and  2.2  summarize  Henrici 
and  Hesse’s  results.  Figure  2.1  shows  that  oil  generally 
lowers  the  heat-transfer  coefficient,  and  that  at  high  heat 
fluxes  and  high  oil  concentrations  (10  percent),  the  effect 
grows  more  pronounced  (slope  decreases).  Figure  2.2  shows 
that  at  some  oil/heat  flux  combinations,  the  heat-transfer 
coefficient  may  actually  be  improved  by  the  addition  of  oil. 
Chongrungreong  and  Sauer  [Bef.  10]  suggest  that  the  heat- 
transfer  behavior  of  refrigerant-oil  mixtures  can  be  attrib¬ 
uted  to  5  major  factors:  1)  the  physical  properties  of  the 

refrigerant-oil  mixture,  2)  the  saturation  temperature  (or 
boiling  pressure) ,  3)  the  tube  diameter,  4)  the  surface 

condition  ox  the  tube  (roughness) ,  and  5)  the  hydrostatic 
liquid  head  above  the  tube. 

1 ♦  Physical  Properties 

Refrigerant- oil  mixtures  have  significantly 
different  physical  properties  than  pure  refrigerants. 
Jensen  and  Jackman  [Bef.  11]  report  that  density  and 
specific  heat  behave  ideally  in  refrigerant-oil  mixtures, 
but  that  viscosity  and  surface  tension  do  not.  Ideal 
behavior  of  refrigerant-oil  mixture  density  and  specific 
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heat  does  not  mean  linear  behavior.  The  governing  egua- 
tions  are: 


1  =  C  _  1  -  C 

Pm  Pot  Prt 


(2.1) 


(1  -  C)  c  +  C  c  „ 
pr  t  po i 


where 


subscripts : 


p  =  density 
C  =  oil  concentration 
cD  =  specific  heat 


(2.2) 


1  =  liquid 
oi  =  mixture 
o  =  oil 

r  »  refrigerant 


Jensen  and  Jackman  report  that  current  refrigerant-oil 
mixture  viscosity  equations  substantially  underpredicted 
their  experimental  data.  No  predictive  equation  has  been 
suggested  for  the  surface  tension  of  refrigerant-oil 
mixtures,  though  Jensen  and  Jackman  developed  a  correlation 
for  R-1 13-oil  mixtures. 

Henrici  and  Hesse  [Ref.  6]  experimentally  determined 
the  surface  tension  for  the  R-1 14-oil  mixtures  that  they 
used  in  their  1971  experiment.  As  shown  in  Figure  2.3,  the 
surface  tension  of  the  mixture  first  decreased  up  to  an  oil 
concentration  of  2.5  percent,  and  then  increased  continu¬ 
ously  with  increasing  oil  concentration.  This  type  of  non¬ 
linear  behavior  makes  explaining  the  change  in  heat- transfer 
coefficient  of  ref rigeraut-oil  mixtures,  due  to  the  changing 
physical  properties  of  these  mixtures,  both  difficult,  and 
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possibly  non-general.  The  behavior  of  refrigerant- oil 
mixtures  is  specific  to  the  particular  mixture  components, 
and  may  be  dependent  upon  the  hind  of  oil  being  used.  Some 
qualitative  consequences  of  adding  oil  to  a  refrigerant, 
however,  can  be  noted. 

The  most  observable  result  of  adding  oil  to  refrige¬ 
rants  is  foaming.  Oil  concentrations  above  1  percent  result 
in  significant  amounts  of  foaming  from  nucleate  boiling. 
The  foam  bubbles  form  because  the  R-114  in  the  R-1 14-oil 
mixture  is  more  volatile  than  oil  and  vaporizes  first, 
creating  a  gas  bubble  surrounded  by  an  oil-rich  layer  (see 
Figure  2.4).  Since  the  bubbles  are  coated  with  oil  film 
with  a  higher  surface  tension  than  the  bulk  liquid  and  a 
have  lower  density,  they  rise  to  the  top  of  the  liquid. 
Because  of  their  surface  tension,  the  bubbles  build  up  on 
the  liquid  surface  to  produce  a  foam  layer. 

This  foaming  action,  which  is  most  pronounced 
between  1-10  percent  oil  concentration  [Ref.  6],  may  affect 
the  heat  transfer  of  tube  bundles  significantly.  For  single 
tubes,  it  is  the  oil  concentration  gradient  which  would  seem 
to  play  the  major  role,  since  the  foam  rises  away  from  the 
tube  surface  and  could  only  interact  with  the  tube  as  it 
sweeps  by  it  from  the  bottom  to  the  top  of  the  tube. 

The  general  decrease  in  the  heat-transfer  coeffi¬ 
cient  upon  adding  oil  to  pure  refrigerants  (recall  Figure 
2.1)  is  subject  to  many  explanations.  Thome  [Ref.  14],  in 
an  extensive  review  of  the  literature,  reports  that  the 
first  explanation  for  the  decrease  in  the  heat- transfer 
coefficient  of  mixtures  was  presented  by  Van  Wijk  et  al.  in 
1956.  The  effect  was  explained  as  being  the  result  of  the 
evaporation  of  the  more-volatile  components,  leaving  an 
oil-rich  layer  with  a  higher  local  boiling  point,  which 
increases  the  amount  of  superheat  required  to  continue 
vaporization  and  bubble  growth,  thus  reducing  the 


Idealized  modal  of  bubble  growth  in  refrigerant-oil  mliture 


Figure  2.4  Oil  Concentration  Gradient  in  a  Bubble 
in  Refrigerant-oil  Mixtures  (from  Ref.  1 t) . 


heat -transfer  coefficient.  Stephan  and  Preusser  [fief.  12] 
demonstrated  conclusively  that  the  work  of  formation  of 
bubbles  in  a  mixture  is  greater  than  in  an  equivalent  pure 
fluid.  They  concluded  that  the  mixture  heat-transfer  coef¬ 
ficient  is  lower  than  for  the  equivalent  pure  fluid  because 
of  the  resulting  decrease  in  the  bubble  population.  Any 
variation  from  a  decreasing  mixture  heat-transfer  coefficent 
(like  Henrici  and  Hesse  show  in  Figure  2.2)  is  attributed  by 
Stephan  [Bef.  13]  to  the  non-linear  variation  of  the  phys¬ 
ical  properties.  Stephan  proposes  that  the  plot  of  surface 
tension  (Figure  2.3),  along  with  thermal-prop erty  varia¬ 
tions,  accounts  for  the  anomalous  rise  in  the  heat- transfer 
coefficient  between  3-6  percent  oil  concentration. 
Chongrongeong  and  Sauer  £Hef.  10]  state  that  it  is  the  rate 
of  heat  diffusion,  governed  by  the  thermal  properties  of  the 
oil-rich  layer,  that  limits  the  bubble  growth  and  that  the 
surface-tension  effects  are  neglible. 

Thome  [Bef.  14]  proposes  that  all  of  the  above 
factors,  as  well  as  the  viscosity  variation,  are  important 
in  explaining  the  rise  in  the  heat-transfer  coefficient  for 
some  refrigerant-oil  mixture  and  heat-flux  combinations. 

In  summary,  all  researchers  agree  that  the  physical 
and  thermal  properties  of  a  refrigerant-oil  mixture  are 
important  factors  in  explaining  the  heat-transfer  behavior 
of  mixtures. 

2 •  Saturation  Temperature 

It  has  long  been  noted  that  increased  saturation 
temperature  (i.e.,  increased  boiling  pressure)  increases  the 
boiling  heat-transfer  coefficient  of  surfaces  in 
refrigerant-oil  mixtures.  In  1963,  Stephan  [Bef.  9]  found 
that  at  high  oil  concentrations,  the  heat-transfer  coeffi¬ 
cient  of  refrigerant-oil  mixtures  becomes  constant  with 
respect  to  the  saturation  temperature.  Stephan  proposed 


that  this  is  because  the  addition  of  oil  to  a  refrigerant 
introduces  a  large  diffusion  resistance,  and  that  since  the 
velocity  of  diffusion  is  almost  independent  of  temperature, 
so  should  the  heat- transfer  coefficient  become  independent 
of  temperature  at  high  oil  concentration.  Figure  2.5  shows 
Henrici  and  Hesse's  data  on  the  effect  of  oil  and  saturation 
temperature  in  R-114.  With  no  oil,  the  effect  of  raising 
the  saturation  temperature  is  seen  to  be  a  rise  in  the  heat- 
transfer  coefficient.  With  oil,  raising  the  saturation 
temperature  is  seen  to  cause  a  slight  drop  in  the  heat- 
transfer  coefficient.  This  effect  has  not  been  explained 
yet. 

3 •  Tube  Diameter 

Cornwell,  Schuller,  and  Einarsson  [Ref.  15]  found 
that  for  smooth  tube  diameters  from  6  mm  to  30  mm,  the 
nucleate  pool  boiling  heat-transfer  coefficient  in  pure 
refrigerants  falls  with  increasing  diameter.  The  effect  of 
tube  diameter  was  correlated  by: 


Nu  =  C  Re2/3 


(2.3) 


where 


Re 


q  d 

U 


Nu 


h  D 
~k~ 


C 


150  for  refrigerants 


Figure  2.5  Effect  of  Saturation  Temperature 
on  B-114-0il  Mixtures  (froa  Bef.  6). 


Since  equation  (2.3)  depends  on  physical  parameters,  the 
effect  of  oil  could  he  significant.  No  known  measurements 
of  the  effect  of  tube  diameter  in  refrigerant-oil  mixtures 
have,  as  yet,  been  performed. 

•  Surface  Condition 

Several  researchers  [Bef.  16,  17,  18]  have  investi¬ 
gated  the  effects  of  surface  roughness  on  the  heat  transfer 
of  pure  refrigerants.  As  surface  roughness  increases,  the 
heat-transfer  coefficient  was  found  to  increase  due  to 
increased  nucleation.  Nishikawa  [Ref.  18]  reported  the 
effect  of  a  variety  of  surface  roughnesses  in  pure  B-114 
over  a  range  of  pressures  from  0.294  MPa  (42  psi)  to 
2. 94  MPa  (420  psi)  . 

5 .  Hydros ta tic  Effect 

The  liquid  column  above  the  boiling  surface  may 
generate  large  static  pressures  which  will  increase  the 
boiling  point.  For  B-114  at  0  oc  (32  °F),  a  0.3  m  (1  ft) 
liquid  head  will  raise  the  saturation  temperature  about  5  °C 
(9  °F) .  In  large  machines,  this  may  be  a  significant 
effect.  However,  for  small  experimental  apparatuses,  the 
effect  is  negligible. 

B.  HtJCLEATE  BOILIHG  OF  BEFBIGEB  ANT-OIL  MIXTURES  FBOM 

ENHANCED  SURFACES 

Webb  [Bef.  19],  in  an  extensive  review  of  the  evolution 
of  enhanced  surface  geometries,  notes  that  the  ability  of 
roughness  to  improve  nucleate  boiling  performance  has  been 
known  for  over  50  years.  However,  it  was  not  until  1968 
that  the  first  commercial  enchanced  surface  was  patented 
[Bef.  20].  Since  then,  the  number  of  commercial  enhanced 
surfaces  has  dramatically  increased  as  the  understanding  of 


their  design  and  operation  has  grown.  Of  the  many  possible 
methods  for  heat- transfer  enhancement,  two  areas  are 
currently  being  commercially  developed:  1)  fins  and  surfaces 
with  reentrant  cavities,  and  2)  porous  coatings. 

1 .  Fins  and  Surfaces  with  Reentrant  Cavities 

The  ability  of  surface  abrasion,  open  grooves,  and 
fins  to  improve  the  heat-transfer  performance  of  a  smooth 
surface  was  first  studied  in  the  1930's.  The  main  diffi¬ 
culty  with  using  surface  abrasion  to  improve  the  heat- 
transfer  performance  is  that  fouling  of  the  surface 
eventually  returns  the  performance  to  that  of  a  non-abraded 
surface.  Studies  in  the  1950's  and  196Q's  centered  on  fins. 
Recent  comparative  studies  f cr  refrigerants  by  Carnavos 
[Ref.  4]  and  Yilmaz  and  Westwater  [Ref.  2]  found  that  fins 
and  grooves  result  in  a  50-100  percent  permament  improvement 
in  the  heat-transfer  performance  compared  to  a  smooth  plain 
tube  in  the  same  refrigerant.  Webb  [Ref.  19],  in  his  liter¬ 
ature  review,  describes  how  researchers  in  the  early  1S60's 
found  methods  to  improve  the  performance  of  fins  by  creating 
reentrant  cavities  on  their  surfaces.  Reentrant  cavities, 
such  as  shown  in  Figure  2.6,  act  as  very  stable  nucleation 
sites  and  thereby  enhance  the  heat-transfer  performance. 
For  a  cavity  to  function  as  a  nucleation  site  and  remain 
active,  even  after  the  surface  is  subcooled,  the  mouth  diam¬ 
eter  (D)  must  fall  within  a  critical  range.  Also,  the 
cavity  must  have  a  reentrant  shape  with  a  maximum  reentrant 
angle  (6).  The  optimum  mouth  diameter  (D)  and  reentrant 
angle  (©)  are  functions  of  the  fluid  properties. 

The  Gewa-T  surface,  patented  in  1979  (manufactured 
by  the  Rieland  Company),  and  the  Thermoexcel-E  surface, 
patented  in  1980  (manufactured  by  the  Hitachi  Company),  are 
two  commercial  surfaces  which  use  modified  fin  shapes  to 
form  the  necessary  reentrant  cavities.  Figure  2.7  shows  the 


(a)  Schematic  cross  section  of  the  C.ewa-T  surface 


(I*)  Schematic  view  of  the  Thermoexcel-E  surfac 


Figure  2.7  Surface  Details  of  Geva- 
Ther»oexcel-E  Reentrant  Surfaces. 
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details  of  these  tube  surfaces.  Carnavos  [fief.  4 ]  found  that 
in  pure  fi-11,  the  Gewa-T  surface  outperformed  a  plain  tube 
by  100-200  percent  and  the  Thermoexcel-E  surface  outper¬ 
formed  a  plain  tube  by  300-400  percent.  Curves  1-5  of 
Figure  2.8  show  the  relative  improvement  in  the  boiling 
heat-transfer  performance  of  B-11  that  can  be  achieved  by 
using  mechanically  produced  reentrant  cavities.  Work 
continues  to  optimize  these  types  of  surfaces  for  the 
various  refrigerants  in  use  commercially.  Both  tubes  have 
been  tested  in  ref rigerant-oil  mixtures  and  did  not  show  a 
significant  decrease  in  performance  [fief.  1  and  7].  The 
cost  of  these  surfaces  is  not  significantly  higher  than  for 
smooth  plain  tubes,  and  the  performance  improvement  is 
dramatic. 

2 .  Porous  Coatings 

The  second  major  type  of  enhanced  surface  is  the 
porous  boiling  surface.  Webb  [fief.  19]  details  the  various 
production  improvements  and  coating  variations  that  have 
been  made  to  the  original  1968  patent  by  Hilton  of  Onion 
Carbide.  The  key  to  the  performance  of  the  porous  coatings 
is  their  small  reentrant  cavities,  which  are  interconnected 
by  substrate  tunnels.  The  particles  used  to  make  the  coat¬ 
ings  are  usually  copper  or  aluminum.  According  to  Webb, 
researchers  have  found  that  the  critical  variable  is  the 
pore  size  rather  than  the  particle  size.  Large  pores  are 
required  for  fluids  with  high  surface  tension  and  high 
thermal  conductivity.  Small  pores  are  optimum  for  fluids 
with  low  surface  tension  and  low  thermal  conductivity  (like 
refrigerants).  Curve  6  of  Figure  2.8  shows  the  relative 
performance  of  the  High  Flux  surface  to  finned  tubes  and 
mechanically  produced  reentrant  surfaces.  Carnavos  [fief.  4] 
found  the  High  Flux  surface  to  be  700-800  percent  better 
than  a  smooth  tube  in  R-11.  No  known  studies  have  been 


published  on  the  performance  of  porous  coatings  in 
refrigerant-oil  mixtures.  Some  studies  of  the  heat-transfer 
performance  of  porous  coatings  in  pure  R-114  and 
refrigerant-oil  mixtures  have  been  made  by  Onion  Carbide, 
but  their  results  are  not  found  in  the  open  literature. 


III.  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 


A.  OVERALL  APPARATUS 

An  overall  schematic  of  the  experimental  apparatus  is 
shown  in  Figure  3.1,  and  a  photograph  is  shown  in  Figure 
3.2.  Karasabun  [Ref.  8]  describes  the  design,  construction, 
and  operation  of  the  apparatus  in  detail. 

The  apparatus  consists  of  two  Pyrex-glass  tees.  Liquid 
R-114  is  boiled  in  glass  tee  (1)  and  is  condensed  in  glass 
tee  (2) .  Gravity  drains  the  condensate  from  the  condenser 
back  to  the  boiling  section.  A  water-ethylene-glycol 
mixture  at  -17  °C  (1  °F)  is  pumped  through  the  condenser 
cooling  coil  via  a  computer-controlled  valve  (VC)  to 
condense  the  R-114  vapor.  The  sump  (7)  that  supplies  the 
water-ethylene-glycol  mixture  is  cooled  by  a  1/2-Ton,  R-12 
air-conditioning  plant. 

Valve  VC  controls  the  R-114  liquid  temperature  and  pres¬ 
sure.  Figure  3.3  is  a  photograph  of  valve  VC  and  the 
computer-controlled  motor  that  operated  VC.  Opening  VC 
causes  more  R-1 14  to  condense  and  lowers  the  system  pres¬ 
sure.  Also,  it  returns  more  subcooled  liquid  to  the  boiling 
section  which  lowers  the  bulk  liquid  temperature.  Since 
data  at  many  heat  fluxes  was  desired  for  a  constant  tempera¬ 
ture,  it  can  be  seen  that  changing  the  heat  flux  without 
adjusting  VC  would  change  the  system  pressure  and  tempera¬ 
ture.  A  computer-controlled  valve  was  thought  to  be  the 
best  way  to  rapidly  return  the  system  to  the  desired  satura¬ 
tion  temperature  following  a  heat  flux  change.  Sections 
III.E  and  IV.  D  describe  in  more  detail  the  computer- 
controlled  valve  and  the  operation  of  the  system  with  the 
computer-controlled  valve  in  use. 
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Oil  was  added  to  the  liquid  R-114  by  draining  it  from  a 
glass  oil  cylinder  (5)  .  The  oil  cylinder  was  refilled  as 
needed  through  valve  7-2  J  m  the  oil  reservoir  (4). 

Two  configurations  of  boiling  tubes  were  tested.  Short 
boiling  tubes  were  tested  to  determine  the  correct  assembly 
procedure  to  obtain  data  on  the  normal  431.8  mm  (17  in.) 
long  boiling  tubes.  The  short  tubes  were  cheaper  to  make, 
thus  more  debugging  attempts  could  be  made  by  testing  them. 
Section  III.C  describes  the  details  of  the  construction  of 
the  short  and  long  test  tubes. 

B.  CI1  SAHPLIHG  APPARATUS 

Following  all  data  runs,  an  attempt  was  made  to  sample 
the  local  oil  concentration  in  the  vicinity  of  a  boiling 
tube.  Figure  3.4  shows  the  oil  sampling  apparatus.  By 
opening  valves  S-1,  S-2,  and  S-4,  the  probe  line  could  be 
purged,  trapping  a  sample  inside  a  flexible  silicon  tube 
that  was  30.5  mm  (12  in.)  long  with  a  3.16  mm  (1/8  in.) 
inside  diameter  using  pinch  clamps.  By  weighing  the  sample 
tube  and  then  boiling  off  the  R-114  leaving  behind  the  oil, 
the  mass  percent  of  oil  in  the  R-1 14-oil  mixture  was  deter¬ 
mined.  The  mass  fraction  of  oil  was  calculated  by: 


Mass  Fraction  =  _ '  (3.1) 

m2  —  ml 

where 

ml  =  mass  of  sample  line 

m2  =  mass  of  sample  line  ♦  R-114  ♦  oil 

m3  =  mass  of  sample  line  ♦  oil  (after  boiling  off  R-114) 


A  Precxsa  Model  80  electronic  mass  balance  was  used  to  weigh 
the  samples.  The  Precisa  Model  80  is  accurate  to 
±  0.  0001  g. 

C.  BGI1IHG  TOBE  COHSTHUCTION 
1 .  Short  Tubes 

Figure  3.5  shows  the  design  of  the  short  tubes.  The 
short  tubes  were  15.9  mm  (5/8  in.)  in  outer  diameter, 
12.7  mm  (1/2  in.)  in  inside  diameter,  and  203.  2  mm  (8  in.) 
in  length.  The  short  tubes  extended  152.4  mo  (6  in.)  into 
the  liquid  R-114  from  the  left  end  flange.  A  25.4  mm  (1  in.) 
long  epoxy  plug  insulated  the  right  end  of  the  tube.  A  1  mm 
(0.04  in.)  thick  copper  disk  behind  the  epoxy  plug  was  soft 
soldered  in  place  to  act  as  a  pressure  barrier.  The  short 
tubes  were  heated  by  a  500-Satt  240-Volt  stainless-steel 
cartridge  heater.  The  heater  was  6.35  mm  (1/4  in.)  in  cuter 
diameter  and  101.6  mm  (4  in.)  in  length. 

The  first  short  tube  was  made  from  thick-walled 
copper  tubing.  This  tube  was  solid  oxygen- free,  high 
conductivity  (OFHC)  copper.  Four  1-2  mm  (3/64  in.)  diameter 
holes  were  drilled  into  the  wall  of  this  tube  at  a  diameter 
of  12.7  mm  (1/2  in.)  for  thermocouple  channels.  Since  this 
tube  was  solid,  and  had  no  sleeve  interface,  it  did  not  have 
an  interface  resistance.  Consequently,  it  was  the  reference 
tube  against  which  all  other  tubes  were  compared  to  deter¬ 
mine  the  amount  of  contact  resistance  they  had. 

Six  other  short  tubes  were  made.  Five  were  made  of 
soft  copper  tubing  and  had  sleeves  inserted  into  the  tube  as 
indicated  in  Table  1.  Soft  soldering  of  the  sleeves  to  the 
tubes  was  determined  to  yield  negligible  contact  resistance 
by  comparison  with  the  solid  tube.  The  last  short  tube  (7) 
was  made  of  90-10  ccpper-nickel  and  was  coated  with  High 
Flux  over  the  active  101.6  mm  (4  in.)  long  section.  This 


Figure  3.4  oil  Sampling  Apparatus. 
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(a)  Sectional  view  of  short  tube. 
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(b)  Thermocouple  sleeve  unwrapped  (at  section  A-A)  to  show 
the  relative  locations  of  the  thermocouple  channels 
(all  dimensions  in  millimeters). 


(c)  Left-end  view  of  the  short  tube. 


Figure  3.5  Sectional  Tiews  of  Short  Boiling  Tube. 


tube  tested  the  compatibility  of  the  sof t-solder- assembly 
method  with  the  copper-nickel  tube  and  High  Flux  coating. 
Section  V.A  describes  the  results  of  data  taken  on  the  short 
tubes  and  the  selection  of  the  soft-solder  method  of 
assembly  for  the  long  tubes. 


TABLE  1 

Summary  of  Short  Tube  Assembly  Hethods 


Tube  Surface 


Eemarks 


1 

2 

3 

4 

5 

6 
7 


Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
High  Flux 


Solid,  thick-walled 
Slide-fit  (0.005  in. 
Slide-fit  (0.002  in. 
Press-fit  0.004  in. 
Press-fit  (0.006  in. 
Sof t-soldered 
Soft-soldered 


tube  (reference) 
clearance) 
clearance) 
interference) 
interference) 


2 .  long  Tube 

Figure  3.6  shows  the  design  of  the  long  boiling 
tubes.  These  boiling  tubes  were  15.9  mm  (5/8  in.)  in  outer 
diameter,  12.7  mm  (1/2  in.)  in  inside  diamter  and  431.3  mm 
(17  in.)  in  length.  The  center  203.2  mm  (8  in.)  was  the 
active  test  section.  For  the  copper-nickel  tube,  the  center 
section  was  the  only  portion  of  the  tube  that  was  coated 
with  High  Flux.  The  remaining  114.3  mm  (4.5  in.)  on  either 
side  of  the  center  section  were  smooth  and  unheated,  and  did 
not  nucleate  under  any  heat  flux  or  oil  condition.  Karasabun 
[Bef.  8]  describes  how  these  end-surfaces  were  treated  by 
the  data-reduction  program  as  an  extended  fin  from  the 
center  section  and  how  their  heat  loss  was  accounted  for. 

The  center  section  was  heated  by  a  1000-Watt 
240-Volt  stainless-steel  cartidge  heater.  The  heater  was 
6.35  mm  (1/4  in.)  in  outer  diameter  and  203.2  mm  (8  in.)  in 
length.  The  heater  was  surrounded  by  a  copper  sleeve  with 
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Figure  3.6  Sectional  ¥iew  of  Long  Boiling  Tube 


eight  1.3  mm  by  1.3  mm  (0.050  in.  by  0.  050  in.)  thermo¬ 
couple  channels  in  them.  Figure  3.7  shows  the  details  of 
the  channel  layout.  The  thermocouple  hot  junctions  were 
welded  to  the  sleeve.  Appendix  A  describes  the  calibration 
of  the  thermocouples.  The  channels  were  oriented  to  provide 
both  axial  and  circumferential  readings  of  the  tube  inner 
wall  temperature.  The  sleevr  was  soft  soldered  to  the  tube. 
The  data  on  the  short  soft-soldered  tube  closely  matched  the 
reference  solid- tube  data,  and  the  short  sof t-soldered-tufce 
data  matched  similiarly  with  long  tube  data.  The  maximum 
circumferential  wall  temperature  variation  in  the  long 
smooth  tubes  was  0.80  K  (0.31  OF)  at  50  kM/m2  compared  to  a 
solid  tube  circumferential  variation  of  0.34  X  (0.61  °F) . 
Section  V.A  describes  in  more  detail  the  circumferential 
variation  of  temperature  that  resulted  using  the  various 
tube  construction  methods.  Some  axial  temperature  variation 
was  experienced  in  the  long  tubes,  particularly  the  long 
High  Flux  tube.  Non-uniform  heat  generation  from  the 
cartridge  heater  is  believed  to  be  responsible  for  the  axial 
temperature  variation.  Section  V.B  descibes  in  more  detail 
the  long  tube  axial  temperature  variation. 

0.  DATA  ACQUISITION  AND  REDUCTION 

A  Hewlett-Packard  3497A  automatic  data  acguistion/ 
control  unit  was  used  to  read  thermocouple  outputs  and  to 
read  an  analog  signal  representing  the  current  and  voltage 
supplied  to  the  cartridge  heater.  A  Hewlett-Packard  9826A 
computer  unit  was  used  to  control  the  HP-3497A  and  to 
analyze  and  store  the  data. 

Information  was  entered  through  the  computer  keyboard  to 
prompt  the  HP-3497A  to  automatically  scan  each  channel.  All 
thermocouple  measurements  were  accomplished  by  0.245  mm 
diameter  (30  gage)  copper-constantan  (type-T)  thermocouples. 
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(b)  Thermocouple  sleeve  unwrapped  (at  section  A-A)  to  show 
the  relative  locations  of  the  thermocouple  channels 
(all  dimensions  in  millimeters). 


LEFT-END  VIEW  RIGHT-END  VIEW 


3  7 


(c)  End  views  of  the  boiling  tube. 

Figure  3.7  Long  Tube  Thermocouple  Channels 


A  power  sensing  device,  described  by  Karasabun  [fief.  8], 
converted  the  AC  current  and  voltage  values  supplied  to  the 
cartridge  heater  to  a  0-5  V  analog,  DC  signal  for  scanning 
by  the  HP-3497  A.  Table  2  lists  the  channel  allocations  in 
the  HP-3497A. 

Following  data  acquisition  for  each  point,  results  were 
computed  according  to  the  step-wise  procedure  outlined  by 
Karasabun  [fief.  8],  and  summarized  in  Appendix  B.  Appendix 
B  also  includes  a  complete  listing  of  the  data-reduction 
program  (DEP2)  . 


TABLE  2 

HP-3 497A  Channel  Allocations 


Channel  Purpose 

25-32  Tube  wall  thermocouples  (location  T(1) 

to  T  (8)  in  long  tube,  T(5)-T(8)  not 
used  in  short  rubes) 

33  Liquid  fi-114  thermocouple  T  (9) 

34  Backup  liquid  K-114  thermocouple  T ( 1 0) 

35  fi-114  vapor  thermocouple  Tplf 

36  Sump  thermocouple  T(12) 

62  Cartridge  heater  voltage  analog  signal 

63  Cartridge  heater  current  analog  signal 


E.  CC  MPUTEB-COHTHOLLED  T1LVE 

The  computer-controlled  valve  (VC)  was  a  Whitney, 
screwed  bonnet,  regulating  valve  with  3/8  in.  Swagelok 
fittings.  The  valve  travel  was  10.5  turns  from  full  shut  to 
full  open.  The  valve  handle  was  replaced  by  a  101.6  mm 
(4  in.)  hard  rubber  disk  (seen  in  Figure  3.3)  which  was 
rotated  by  the  motor  pinion  gear. 

The  computer-controlled  motor  was  a  General  Electric 
"Minigear  Motor”  with  a  speed  of  105.7  rpm  and  a  torque  of 
3.39  N-m  (30  lb-in) .  The  motor  direction  was  controlled  by 
two  sets  of  Crydom  solid-state  relays  that  acted  to  open. 
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shut,  or  hold  the  motor  depending  on  signals  from  the 
computer.  The  computer  controlled  the  amount  of  time  that 
the  valve  vas  moving.  Approximately  3  seconds  were  required 
for  the  motor  to  turn  the  valve  one  turn. 

The  flow  through  the  valve  vas  checked  by  a  flow  meter 
and  verified  to  be  approximately  linear  (1  turn  =  10  percent 
flow)  for  the  flows  most  often  used.  The  valve  position  was 
tracked  by  a  10- turn  potentiometer,  which  was  connected  to 
the  computer-controlled  motor  by  another  101.6  mm  (4  in.) 
hard  rubber  disk  driven  off  the  motor  pinion  gear.  Hard 
rubber  disks  were  used  instead  of  metal  gears  to  avoid 
damage  to  the  gear  teeth  during  the  program  debugging  stage. 
The  rubber  disks  allowed  sufficient  slipping  when,  for 
example,  the  computer  sent  a  valve-open  signal  even  though 
the  valve  was  fully  cpen. 

The  B-114  liquid  temperature  input  to  the  computer  was 
provided  by  a  separate  copper-constantan  thermocouple 
installed  in  the  same  liquid  B-114  thermocouple  well  that 
the  HP-3497A  data  acquistion/control  unit  used.  The  B-114 
theromocouple  emf  for  the  computer  was  amplified  by  an  Omega 
thermocouple  DC  millivolt  amplifier  before  input  in  to  the 
computer.  The  data  acquisition  system  and  the  computer- 
controlled  valve  system  were  completely  independent  systems. 

The  computer  used  for  the  computer-controlled  valve  was 
an  Octagon  Systems  SYS-2A  microcomputer  with  an  Esprit  I 
terminal  connected  via  an  HS-232C  serial  port.  Appendix  C 
lists  the  control  program  used.  The  control  program  vas 
written  in  NSC  "Tiny  EASIC.  "  The  control  algorithm  simulates 
a  proportional-integral-derivative  (PID)  controller  to  vary 
the  valve  opening  and  shutting  times.  In  the  final  program, 
limits  were  placed  on  the  numerical  value  of  certain  program 
variables  to  prevent  register  overflow,  jamming  the  valve 
fully  open/shut,  and  to  lessen  the  impact  of  system  noise  on 
the  response  of  the  system.  NSC  Tiny  BASIC  is  limited  to 


integers  from  -32000  to  32000  and  fractional  numbers  are 
truncated.  The  values  of  the  weighting  factors  (A,B,C)  for 
the  proportional  ten  (E)  ,  integral  term  (I)  ,  and  derivative 
term  (D)  were  determined  by  trial-and-error. 

The  control  algorithm  simplifies  to  the  following  lines: 


60  Input  required  temperature  E 
230  Read  E-1T4  temperature  M 


(D) 


400  Valve  command  V=  (E/A)  ♦  (B^D)  +  (I/C) 
420  If  (V>0)  then  open  valve 
430  Otherwise  shut  valve 
540  GOTO  230 


error  is  positive 
f  error  is  negative 


Section  V.C  desrites  in  detail  the  system  response  under 
this  algorithm.  Computer  control  of  valve  VC  was  not  used 
during  the  data  taking  as  originally  planned.  Section  V.C 
discusses  the  reasons  for  this  decision. 


17.  E1P  BRIBE  IT  AI  PBOCBDOBZS 


A.  I1STA1.LATIOI  Of  TOBE  II  APPARATUS 

Prior  to  installation  in  the  apparatus,  the  boiling  tube 
surface  was  cleaned  with  a  2  percent  Nital  solution  (to 
remove  surface  oxidation  and  oil),  rinsed  with  acetone,  and 
air  dryed.  Short  tubes  tested  with  and  without  the  above 
treatment  showed  no  change  in  the  heat-transfer  performance. 
The  treatment  was  effective  in  removing  the  slight  surface 
oxidation  present  following  soft  soldering  without  changing 
the  smooth  or  High  Flux  tube  performance.  Additionally, 
following  testing  with  oil,  a  treatment  was  needed  to  return 
the  High  Flux  surface  to  the  "no-oil"  condition  for  further 
testing. 

After  installing  the  tube  in  the  glass  tee,  the  appa¬ 
ratus  was  evacuated  to  29  in.  Hg  by  the  portable  mechanical 
vacuum  pump  (6)  shown  in  Figure  3.1.  System  pressure  was 
measured  by  a  Marsh  pressure  gage  (30  in.  Hg  to  150  psi 
range,  ±2.5  in-  Hg  and  ±0.5  psi  accuracy).  The  apparatus 
was  left  at  vacuum  for  2  hours  to  check  for  leaks  prior  to 
each  run.  No  noticabl^  vacuum  was  observed  within 
the  accuracy  of  the  .  Next,  the  system  gage 
pressure  was  rais  Mr-  27  psi),  the  saturation 
pressure  of  8-114  ,  by  opening  valve  77  (see 
Figure  3.1  for  t*  -  c  on  of  the  valves).  An 
Automatic  Halogen  Leav  _,  TIF  5000,  was  used  to  check 
for  H-114  leakage.  s  -sitivity  of  this  detector  is  3 
ppm  minimum  concentr/  on.  After  pressure  egualization  with 
the  P.-114  reservoir  ^3)  ,  the  reservoir  drain  valve  (76)  and 
condenser  return  valve  (75)  were  opened  to  fill  the  boiling 
tee  with  liquid  fi-114. 


Prior  to  installation,  the  left  end  flange  had  been 
scribe  marked  to  indicate  the  liguid  level  corresponding  to 
1600  cc  (2500  gm)  of  liguid  R-114  at  21  °C  (70  °F)  .  This 
was  the  mass  of  pure  R-114  in  the  apparatus  at  the  beginning 
of  the  data  runs.  The  apparatus  was  now  ready  for  taking 
data . 

B.  GENERAL  OPERATION 

Table  3  lists  the  109  data  runs  accomplished  during  this 
thesis  effort  and  their  purpose.  The  data  runs  were 
numbered  seguentially  and  preceeded  by  a  2-  or  3-letter 
prefix  to  indicate  the  tube  type.  The  tube  prefixes  were: 

SS  =  Short  Solid  Tube 
HF  =  High  Flux  Long  Tube 

WH  =  Wieland  Hard  Copper  Long  Smooth  Tube 

SSF  =  Short  Slide- Fit  Tube 

SPF  =  Short  Press-Fit  Tube 

SST  =  Short  Soft-Soldered  Tube 

SHF  =  Short  High  Flux  Tube 

The  short  tube  runs  consisted  of  6  data  points  at  6 
different  heat  fluxes  (usually  59,  37,  22,  14,  8,  and  5 
kN/m2).  The  normal  tube  runs  consisted  of  at  least  7 
different  heat  fluxes  with  6  data  readings  at  each  heat  flux 
(usually  98,  6  1,  37,  22,  14,  8,  and  5  kK/m2)  with  some  addi¬ 
tional  low  heat  fluxes  investigated  for  the  0,  3,  and  10 

percent  oil  cases  to  check  for  the  onset  of  nucleate  boiling 
and  hysteresis  in  the  High  Flux  surface. 

In  all  cases,  the  data  set  was  begun  by  starting  the 
cooling  pump  (8)  and  opening  valve  VC  slightly  to  slowly 
cool  the  liguid  R-114  to  the  desired  saturation  temperature. 
The  R-114  vapor  temperature  was  also  monitored  and  when  both 
liguid  and  vapor  temperatures  were  stable  (usually  after 
30-45  minutes)  ,  the  heat  flux  would  be  established,  the 
saturation  temperature  reestablished,  and  the  data  aguisi- 
tion  unit  was  allowed  to  take  data.  Usually,  the  vapor 
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TABLE  3 

Summary  of  Data  Bans 


Run 

Tsat 

No. 

No. 

(°C) 

Pts 

WH10 

SPF  1 1 

6.7 

3 

SPF12 

-2.2 

6 

SPF  13 

2.2 

6 

SPF14. 

6.7 

6 

SSF15 

-2.2 

6 

SSF16 

2.2 

6 

SSF17 

6.7 

6 

SSF18 

-2.2 

6 

SSF19 

-2.2 

6 

SSF2  0 

-2.2 

6 

SSF21 

-2.2 

6 

SSF22 

6.7 

6 

SSF23 

-2.2 

6 

SPF  2  4 

-2.2 

6 

SPF25 

-2.2 

6 

SPF26 

-2.2 

6 

SPF27 

-2.2 

6 

SPF  2  8 

6.7 

6 

SSF29 

-2.2 

6 

SSF30 

6.7 

6 

SPF31 

-2.2 

6 

SPF3  2 

-2.2 

6 

SPF33 

2.2 

6 

SPF34 

6.7 

6 

SS35 

-2.2 

6 

SS36 

-2.2 

6 

SS37 

6.7 

6 

SPF  3  8 

-2.2 

6 

SPF39 

-2.2 

6 

SPF  40 

-2.2 

6 

SPF  41 

-2.2 

6 

SS42 

-2.2 

6 

SST43 

-2.2 

6 

SST44 

-2.2 

6 

SHF45 

-2.2 

4 

SHF46 

-2.2 

5 

SHF  4  7 

-2.2 

7 

SHF48 

-2.2 

8 

SHF  4  9 

-2.2 

6 

SHF50 

-2.2 

6 

SHF  51 

-2.2 

6 

SHF52 

-2.2 

6 

SHF53 

-2.2 

6 

SHF  54 

-2.2 

6 

SHF55 

-2.2 

6 

SHF56 

-2.2 

6 

SHF57 

-2.2 

6 

SHF  5  8 

-2.2 

6 

SHF59 

-2.2 

6 

SHF60 

-2.2 

6 

SHF6  1 

-2.2 

6 

SHF  6  2 

-2.2 

6 

SHF63 

-2.2 

6 

SHF64 

-2.2 

6 

SHF65 

-2.2 

6 

SHF  6  6 

-2.2 

6 

SHF67 

-2.2 

7 

SHF68 

-2.2 

6 

Remarks 

Buns  JJH01  to  NH 10  taken  by  Karasabun 
Debug  nev  tube 

Press-fit  tube  4  (0.004  in.  interference) 
Effect  of  Tsat 
Effect  of  Tsat 

Slide-fit  tube  2  (0.005  in.  clearance) 

Effect  of  Tsat 

Effect  of  Tsat 

Rotate  tube  90  degrees 

Study  low  beat  flux  error  band 

Study  high  beat  flux  error  band 

Rotate  heater  -90  degrees,  tube  fixed 

Shift  heater  and  thermocouples 

Repeatability 

Repeatability 

Repeatability 

Shift  heater  and  thermocouples 

Clean  with  Nital,  acetone,  and  air  dry 

Repeatability 

Slide-fit  tube  3  (0.002  in.  clearance) 

Effect  of  Tsat 

Repeatability 

Time  Variation  Study  (1  day  later) 

Effect  of  Tsat 
Effect  of  Tsat 
Short  Solid  Tube 
Repeatability 
Effect  of  Tsat 
4  days  later 

Shift  heater  and  thermocouples 
7  days  later 
10  days  later 

Time  Study  Solid  Tube  (no  effect) 

Short  soft- soldered  tube  6 
Repeatability 

Debug  short  High  Flux  tube  7 
Study  Tvapor  superheat 
Repeatability 
Repeatability 

Study  hydostatic  head  (+1  in.  level) 
-Increase  data  pts  at  heat  fluxes 
-Runs  49-54  form  data  set 
-Each  run  at  different  heat  flux 
-Compare  with  set  55-60 
End  of  Set  49-54 

Study  hydrostatic  head  (+0.5  in.  level) 

-same  as  above  set  49-54 

-same 

-same 

-same 

End  of  Set  55-60 

Install  Tvapor  radiation  shield 

Study  Tliquid  subcooling 

Clean  High  Flux  w/Nital  and  acetone 

Repeatability 

Repeatability 

Repeatability 

Rotate  tube  +90  and  *180  degrees 
Repeatability 


Table  3 

Susaary  of  Data  Runs  (coni' d) 


Run 

Tsat 

No. 

No. 

(°C) 

Pts 

SHF  6  9 

-2.2 

7 

SHF70 

6.7 

36 

SHF71 

-2.  2 

36 

SHF72 

-2.2 

36 

SHF73 

-2.2 

36 

SHF74 

-2.2 

36 

SHF75 

-2.2 

36 

SHF76 

-2.2 

36 

SHF77 

-2.2 

36 

RH78 

-2.2 

9 

HH79 

-2.2 

42 

RH80 

-2.2 

1 1 

WH81 

-2.2 

42 

NH8  2 

-2.  2 

42 

WH83 

-2.2 

38 

RH84 

6.7 

49 

TJH85 

-2.2 

66 

HH86 

-2.2 

48 

WH87 

6.7 

48 

WH88 

-2.2 

48 

HH8  9 

6.7 

48 

WH90 

-2.  2 

36 

WH9 1 

-2.2 

48 

WH92 

6.7 

48 

RH93 

6.7 

48 

WH9  4 

-2.2 

48 

NH95 

-2.2 

61 

HH96 

6.7 

53 

HH97 

-2.2 

66 

RH98 

-2.2 

45 

»H99 

6.7 

66 

HP  100 

-2.2 

66 

HF101 

-2.2 

72 

HF102 

-2.2 

24 

HF103 

6.7 

66 

HF104 

6.7 

69 

HF105 

6.7 

56 

HF106 

-2.2 

48 

HF107 

6.7 

48 

HF108 

6.7 

48 

HF109 

-2.2 

48 

HF1 10 

6.7 

48 

HF111 

-2.2 

39 

HF112 

-2.2 

72 

HF113 

-2.2 

48 

HF1 14 

6.6 

48 

HF115 

-2.2 

48 

HF116 

-2.2 

72 

HF117 

6.7 

72 

HF118 

6.7 

48 

SPF1  19 

-2.2 

36 

Re sarks 

Apply  thermal  grease  to  heater  (no  effect) 

Effect  of  Tsat 

Repeatability 

Debug  oil  addition  (0.2  percent  added) 
Remove  and  clean  tube,  aad  1  percent  oil 

2  percent  oil 

3  percent  oil 
6  percent  oil 

10  percent  oil 

Axial/Circumferential  variation  study 
0  percent  oil.  decreasing  g 
Rotated  tube  -90  and  -180  degrees 
Repeatability  (of  run  HH79) 

0  percent  oil,  increasing  g 
Repeatability  (of  run  WH82) 

0  percent  oil,  increasing  g 
0  percent  oil,  decreasing  g 
1  percent  oil,  decreasing  g 

1  percent  oil,  decreasing  g 

2  percent  oil,  decreasing  g 

2  percent  oil,  decreasing  g 

3  percent  oil,  increasing  g 
3  percent  oil,  decreasing  g 
3  percent  oil,  decreasing  g 
6  percent  oil,  decreasing  g 
6  percent  oil,  decreasing  g 

10  percent  oil,  increasing  g 
10  percent  oil,  increasing  g 
10  percent  oil,  decreasing  g 
Repeatability  (of  run  WH9  5) 

10  percent  oil,  decreasing  g 
0  percent  04. 1,  increasing  g 
0  percent  oil,  decreasing  g 
Repeatability  (of  run  HF100) 

0  percent  oil,  increasing  g 
0  percent  oil,  decreasing  g 
Repeatability  (of  run  HF104) 

1  percent  oil,  decreasing  g 

1  percent  oil,  decreasing  g 

2  percent  oil,  decreasing  g 

2  percent  oil,  decreasing  g 

3  percent  oil,  decreasing  g 
3  percent  oil,  increasing  g 
3  percent  oil,  decreasing  g 
6  percent  oil,  decreasing  g 
6  percent  oil,  decreasing  g 

10  percent  oil,  increasing  g 
10  percent  oil,  decreasing  g 
10  percent  oil,  decreasing  g 
10  percent  oil,  increasing  g 
Complete  study  of  time  variation  SPF 


Rote:  75  data  runs  for  Computer-Controlled  Valve 
testing  not  listed. 


temperature  would  read  higher,  up  to  2  K,  than  the  liquid 
temperature  due  to  the  vapor  becoming  superheated.  The 
apparatus,  though  coded  by  the  B-114,  was  still  hot  enough 
(60  °F)  to  superheat  the  vapor.  Measurements  of  the  liquid 
temperature,  system  pressure,  and  vapor  temperature  (with 
vapor  probe  shielded  by  a  radiation  shield)  confirmed  this. 
The  liquid  temperature  best  indicated  the  saturation 
temperature. 

The  HP-3497A  data  aguisition  unit  would  scan  each 
channel,  compute  the  heat- transfer  coefficient,  print  the 
results  (an  example  printout  is  shown  in  Appendix  D)  ,  and 
store  the  data  on  the  floppy  disk. 

Following  the  taking  of  all  data  points,  the  data  set 
was  statistically  analyzed  by  subroutine  STATS  to  compute 
the  average  heat-transfer  coefficient  at  each  heat  flux  and 
the  standard  deviation  of  the  6  data  points  for  a  given  heat 
flux.  The  standard  deviation  of  the  6  data  points  was 
usually  0.5  percent  for  the  heat  flux  and  1-2  percent  for 
the  heat-transfer  coefficient. 

After  all  data  sets  at  a  given  oil  concentration  were 
complete,  oil  was  added  via  valve  71.  The  oil  immediately 
dissolved  in  the  B-114.  No  carryover  to  the  condenser  was 
noted  except  for  several  small  drops  at  10  percent  oil  and 
the  highest  heat  flux  during  the  last  few  data  sets. 
Foaming  occurred  with  the  addition  of  oil,  and  increased 
with  both  increasing  heat  flux  and  increasing  oil  concentra¬ 
tion.  Figure  4.1  to  4.3  show  photographs  at  heat  fluxes  of 
30  kW/m2  and  98  kw/m2  for  oil  concentrations  of  0,  3,  and  10 
percent.  When  oil  addition  was  not  taking  place,  the  oil 
cylinder  and  reservoir  were  isolated  from  the  apparatus  by 
valves  V3  and  74  to  minimize  B-114  absorption  by  the  oil. 


(a)  3  percent  oil  at  37  kW/in^. 


Pllff 


C.  S AHPLIIG  OIL 


Oil  sampling  was  done  at  room  temperature  to  prevent  the 
temperature  of  the  sampling  lines  from  boiling  off  the  B-114 
and  yielding  false  readings.  The  sampling  procedure  used 
was: 

1.  Weigh  the  empty  flexible  sample  tubing  and  pinch 
clamps.  Connect  them  to  the  sampling  apparatus. 

2.  Open  S- 1 ,  S-4,  and  set  S-2  to  position  2  to  provide  a 
purge  path  for  the  sample  flow. 

3.  lower  the  sample  container  to  a  height  below  the 
level  of  the  B-114  glass  tee. 

4.  After  purging  the  sample  probe  and  lines,  use  pinch 
clamps  to  isolate  the  304.8  mm  (12  in.)  flexible 
sample  tubing  and  trap  a  sample.  Shut  S-1  and  switch 
S-2  to  position  1. 

5.  Disconnect  the  flexible  sample  tubing  and  immerse  it 
in  ice  to  lower  its  internal  pressure  below  atmos¬ 
pheric  pressure. 

6.  Open  S-3  to  provide  a  vent  for  the  purge  container 
and  pour  the  purge  volume  back  into  the  boiling  tee 
by  lifting  it  to  a  height  so  that  it  flows  by 
gravity . 

7.  Weigh  the  sample  line  (the  mass  balance  used  was 
accurate  to  ±  0.0001  g)  . 

8.  Open  a  pinch  clamp  and  allow  air  to  warm  the  sample 
line  and  evaporate  off  the  liguid  E-114.  Keep  the 
sample  line  on  the  mass  balance  to  collect  any  drops 
of  oil  that  may  splatter  as  the  liguid  B-114  evapo¬ 
rates. 

9.  After  allowing  several  hours  for  the  R-114  to  evapo¬ 
rate,  reweigh  the  flexible  sample  tubing. 

10.  Calculate  the  mass  percent  of  oil  using  equation 


V.  SYSTEM  OPERATION  AND  PROBLEMS 


A.  CIRCUMFERENTIAL  NALL  TEMPERATURE  VARIATION 

Following  the  construction  of  the  solid,  thick- walled, 
reference  tube,  an  investigation  was  begun  to  develop  an 
assembly  method  which  would  result  in  negligible  contact 
resistance  between  the  inner  tube  wall  and  the  copper 
sleeve,  and  yield  data  comparable  to  that  for  the  solid 
tube. 

Short  slide-fit  tubes  2  and  3  (see  Table  1)  had  sleeves 
that  could  be  slid  into  the  test  tube  with  tube  3  having  a 
tighter  clearance.  The  slide-fit  tubes  exhibited  a  circum¬ 
ferential  wall  temperature  variation  of  7-14  K  at  50  kW/m2 
in  pure  R-114  boiling  at  -2.2  oc.  This  variation  matched 
data  on  a  similar  long  tube  (with  slide-fit  sleeve)  tested 
by  Karasabun  [Hef.  8].  Tube  2  had  an  average  wall  tempera¬ 
ture  of  33  °C  and  tube  3  had  an  average  wall  temperature  of 
25  °C  at  50  kN/m2.  The  tighter  clearance  of  tube  3  resulted 
in  a  lower  contact  resistance,  and  the  wall  temperature 
subseguently  dropped.  The  short  solid  tube  under  similar 
conditions  exhibited  a  0.34  K  circumferential  wall  tempera¬ 
ture  variation  and  an  average  wall  temperature  of  only 
10.7  °C .  Sauer  et  al.  [Ref.  21  with  a  similar  experi¬ 
mental  apparatus,  used  a  mechanically-press-fitted  brass 
sleeve  for  obtaining  data.  Tubes  4  and  5  had  copper  sleeves 
that  were  mechanically  cold  pressed  into  the  tube.  The 
diametral  interference  was  0.01  mm  (0.0004  in.)  for  tube  4 
and  0.015  mm  (0.0006  in.)  for  tube  5.  The  interface  pres¬ 
sure  obtained  was  calculated  to  be  15.2  MPa  (2200  psi)  for 
tube  4  and  22.8  MPa  (3300  psi)  for  tube  5.  During  pressing, 
the  sleeves  were  lubricated  with  glycerin  (C3  ft  O3  )  a 


high-thermal-conductivity  compound  that,  according  to 
Incropera  and  Dewitt  £Bef.  22  ]#  should  result  in  a  10-100 
times  reduction  in  the  contact  resistance  when  combined  with 
the  contact  pressures  mentioned  above. 

Tubes  4  and  5  had  circumferential  wall  temperature  vari¬ 
ations  of  0.9-1. 5  K  with  an  average  wall  temperature  of 
16.2  °C  for  tube  4  and  14.7  °C  for  tube  5  during  initial 
testing.  However,  with  time,  the  circumferential  wall 
temperature  variation  grew  to  5-8  K  and  the  average  wall 
temperature  increased  to  about  28  °C  for  both  tubes.  The 
resulting  drop  in  heat-transfer  performance  is  shown  in 
Figure  5.1.  This  result  is  believed  to  be  due  to  the  phenom¬ 
enon  referred  to  as  "stress  relaxation." 

Stress  relaxation  is  a  form  of  creep.  The  Metals 
Handbook  [Bef.  23]  notes  that  copper  alloys  easily  undergo  a 
decrease  in  stress  resulting  from  transformation  of  elastic 
strain  into  plastic  strain  in  a  constrained  solid.  The 
phenomena  occurs  even  at  relatively  low  operating  tempera¬ 
tures  (at  25  °C ,  an  80  percent  drop  in  stress  can  occur  in 
200  hours) ,  and  is  of  most  concern  in  applications  like 
press-fits  and  solderless  wrapped  copper  connectors.  The 
time- dependent  data  of  Figure  5. 1  appears  to  be  the  result 
of  a  dropping  press-fit  interface  pressure  yielding  a  higher 
contact  resistance,  higher  inner  wall  temperatures,  and 
lower  heat-transfer  coefficients. 

Stephan  and  Mitrovic  £Hef.  7]  used  a  combination  of 
mechanically  press-fit  and  soft-soldered  inner  sleeves  to 
obtain  data  on  the  Gewa-T  surface  in  R-114.  Soft-soldering 
was  initially  overlooked  by  this  experimenter  because  the 
normal  heating  method  used  is  an  oxy-acet ylene  torch  on  the 
tube  surface.  The  flame  temperatures  of  an  oxy-acet ylene 
torch  are  above  800  ®C  (1500  °F)  .  The  resulting  oxidation 
and  heat  damage  (the  High  Flux  coating  melts  at  800  °C)  to 
the  High  Flux  surface  was  unacceptable.  Tests  showed 


however#  that  the  inner  cartridge  heater  could  be  used  to 
apply  a  controlled  amount  of  heat  to  melt  a  low  temperature 
(melting  point  160  og  solder#  while  closely  monitoring  the 
tube  temperature  using  the  sleeve  thermocouples.  This 
method  was  used  to  produce  soft-soldered  short  tube  6  which 
showed  circumferential  wall  temperature  variations  of  0. 8  K 
and  an  average  wall  temperature  of  11.5  °C  at  50  kfl/m2. 
Figure  5.2  compares  the  results  of  the  soft-soldered  and 
solid  tubes.  The  agreement  between  the  tubes  is  excellent. 
The  long  boiling  tubes  matched  the  short  soft-soldered  data 
very  closely. 

Tests  made  by  rotating  the  test  tubes  90  degrees  and  180 
degrees  showed  that  the  slight  circumferential  wall  tempera¬ 
ture  variation  of  these  tubes  was  due  to  the  surface  charac¬ 
teristics  of  the  boiling  tube  rather  than  due  to  the 
thermocouples.  Thermocouples  located  near  more  active 
nucleation  sites  of  the  smooth  and  High  Flux  surfaces  had 
slightly  lower  local  wall  temperatures.  The  0.8  K  circum¬ 
ferential  wall  temperature  variation  that  resulted  from 
soft -soldering  the  sleeves  of  the  tubes  is  much  smaller  than 
the  11  K  variation  reported  by  Sauer  et  al.  [Bef.  21],  and 
is  about  the  same  as  the  1  K  variation  reported  by  Stephan 
and  Mitrovic  [Bef.  7]. 

B.  AIIA1  TEH PE BAT OB E  VARIATION  OF  LONG  TUBES 

Karasabun  [Bef.  8]  reported  an  axial  temperature  varia¬ 
tion  of  about  20  K  along  the  inner  wall  of  his  slide-fitted 
sleeve.  Stephan  and  Mitrovic  [Bef.  7]  reported  an  axial 
temperature  variation  of  1  K  along  the  sleeve  they  used. 
The  long  soft-soldered  tubes  tested  in  this  experiment 
exhibited  an  axial  temperature  distribution  that  varied  with 
heat  flux.  Figures  5.3  and  5.4  show  the  axial  temperature 
distribution  as  a  function  of  position  along  the  boiling 


Tube 


surface  for  the  smooth  tube  and  High  Flux  tube  respectively. 
The  variation  at  20  kW/m2  (the  heat  flux  at  which  Stephan 
and  Mitrovic  cite  their  1  K  variation)  is  less  than  1  K  for 
both  tubes.  The  maximum  axial  temperature  variation  is  3  K 
for  the  High  Flux  surface  at  98  kH/m2. 

The  cartridge  heaters  used  were  precision-wound, 
ceramic-core,  magnesium-oxide-insulated,  Incoloy-sheathed 
WATLOS  FIBEEOD  heaters.  The  heaters  had  6.35  mm  (0.25  in.) 
long  lava-rock  plugs  at  either  end  of  the  heater  that 
reduced  the  actual  heating  length  on  each  end.  The  heaters 
were  initially  believed  to  generate  a  uniform  heat  flux  at 
all  power  settings,  but  the  axial  wall  temperature  data 
indicate  the  heat  generation  varied  with  the  power  level. 

Since  the  axial  temperature  distribution  was  fairly 
linear  over  most  practical  heat  fluxes  (less  than  37  kW/m2)  , 
the  arithmetic  average  of  all  8  wall  thermocouples  was  used 
to  conpute  the  heat- transfer  coefficient.  This  results  in 
slightly  lower  heat- transfer  coefficients  and  is  a  conserva¬ 
tive  estimate  of  the  performance  of  the  High  Flux  surface. 
Appendix  E  analyzes  the  resulting  error  in  the  heat-transfer 
coefficient  from  using  the  arithmetic  average  for  calcu¬ 
lating  the  heat-transfer  coefficient. 

C.  CCHPUTEB-CONTHOL1ED  VALTE 

Figure  5,5  shows  the  system  response  to  closing  valve  VC 
manually  and  the  improved  response  of  the  system  when  using 
the  computer-controlled  valve.  3y  correctly  cycling  the 
valve  between  open  and  shut,  the  system  saturation  tempera¬ 
ture  could  rapidly  be  changed.  Establishing  a  stable  egui- 
librium  temperature  after  a  large  valve  movement  was, 
however,  difficult  during  the  initial  testing  of  the  appa¬ 
ratus.  The  extensive  trial-and-error  testing  for  the 
correct  weighting  factors  of  the  proportional,  integral,  and 
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derivative  terms  resulted  in  a  large  data  base  concerning 
the  saturated  system  operating  characteristics  and  a  large 
amount  of  practice  for  the  operator  in  manual  control  of 
valve  VC  following  unsuccessful  computer-controlled 
transients. 

Manual  operation  of  valve  VC  was  done  in  conjunction 
with  the  section  of  the  data- reduction  program  (lines  1890 
to  2485  of  program  DEF2  listed  in  Appendix  B)  that  monitored 
the  saturation  temperature  and  the  rate  of  change  of  satura¬ 
tion  temperature  continuously.  It  was  found  that,  because 
of  the  better  temperature  resolution  of  the  HP-3497A  data 
acquisition  unit,  with  manual  operation  and  a  trained  oper¬ 
ator,  temperature  transients  could  be  maintained  within 
±  0.02  K.  This  was  a  much  tighter  control  band  than  was 
possible  with  the  computer-controlled  valve  because  of  the 
poor  resolution  of  its  temperature  sensing  input  from  the 
Omega  thermocouple  DC  millivolt  amplifier.  Additionally, 
after  much  practice  and  experience,  the  manual  method  was 
found  to  be  as  fast  as,  or  faster  than  the  use  of  computer- 
controlled  valve. 

To  obtain  the  highest  accuracy  data  in  the  shortest 
period  of  time,  computer-control  of  valve  VC  was  abandoned. 
Section  VII.B  includes  recommendations  for  further  improve¬ 
ments  to  the  computer  control  system  to  restore  its 
usefulness. 
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(a)  Natural  system  response  to  shutting  valve  VC  to  cause 
a  desired  increase  in  system  temperature. 


TIME  CM IN) 

(b)  Improvement  in  system  response  time  using  computer 
control  of  valve  VC  to  reach  a  desired  temperature 
(note  the  different  time  scale). 


Figure  5.5 


Systen  Response  to  a  Transient  with 
Computer-Controlled  Valve. 


71.  BESPITS  AND  DISCUSSION 


A.  BGILIBG  PERFORMANCE  OF  SHOOTO  TOBE 

Figure  6.1  shows  the  nucleate  pool-boiling  performance 
of  the  smooth  copper  tube  in  B-114.  The  behavior  with  no 
oil  represents  typical  nucleate  pool-boiling  performance. 
From  point  A  to  point  B,  a  continuous  increase  in  wall 
superheat  (TWQ-Tsat)  is  observed  when  heat  flux  is 
increased.  No  bubbles  were  observed  in  this  region  of  the 
curve  as  this  region  represents  natural  convection.  At 
point  B,  incipient  nucleate  boiling  occurs.  From  point  3  to 
point  C  (or  C'  for  10  percent  oil),  a  reduction  in  wall 
superheat  is  observed  while  the  heat  flux  is  continuously 
increased.  This  region  is  known  as  the  mixed-boiling 
region,  where  transition  from  natural  convection  to  nucleate 
pool  boiling  takes  place.  In  this  region,  the  heated 
portion  of  the  tube  began  to  activate  an  increasing  number 
of  nucleation  sites,  while  the  unheated  ends  showed  no 
bubbles.  In  fact,  the  unheated  ends  underwent  only  natural 
convection,  due  to  axial  conduction  of  heat  along  the  tube 
wall,  at  all  heat  fluxes.  The  transition  from  natural 
convection  to  nucleate  boiling  occurred  rapidly  when  there 
was  no  oil  present.  The  surface  would  burst  into  nucleate 
boiling  in  less  than  a  second  after  the  first  nucleation 
site  became  active.  At  point  C  (or  C’),  all  the  available 
nucleation  sites  were  apparently  active.  After  point  C,  the 
wall  superheat  again  increases  with  increasing  heat  flux  as 
shown  in  region  C-D.  In  region  C-D,  no  new  nucleation  sites 
were  seen  to  become  active.  Instead,  the  bubble  departure 
rate  increased.  When  the  heat  flux  is  decreased  after  having 
established  complete  nucleate  boiling,  the  curve  follows  a 


different  path  (point  D  to  point  E)  .  The  existence  of 
stable  nucleation  sites,  which  remain  active  over  a  wide 
range  of  heat  fluxes,  results  in  better  heat- transfer 
performance  than  in  natural  convection,  resulting  in  a  lower 
wall  superheat. 

The  effect  of  adding  oil  is,  according  to  Stephan 
[Ref.  9]#  to  introduce  a  mass  diffusion  resistance  and  lower 
the  heat-transf er  coefficient.  As  seen  in  Figure  6.1,  in 
region  A-B,  both  the  3  and  10  percent  curves  were  lower  than 
the  0  percent  oil  curve.  The  3  percent  oil  curve  is  lower 
than  the  10  percent  oil  curve  probably  because  of  the  non¬ 
linear  physical  property  variations  of  refrigerant-oil 
mixtures.  The  non-linear  variation  of  surface  tension  (see 
Figure  2.3}  would  not  seem  to  be  responsible  for  this 
anomally.  The  curves  in  region  A-B  support  the  contentions 
of  Chongrungeong  and  Sauer  [Ref-  10]  and  Thome  [Bef.  14] 
that  the  non-linear  variation  of  physical  properties  of 
refrigerant-oil  mixtures,  other  than  surface  tension, 
explains  the  heat-transfer  behavior  of  refrigerant-oil 
mixtures.  The  effect  of  adding  10  percent  oil  was  to  delay 
the  transition  to  complete  nucleate  boiling  on  the  tube. 
With  10  percent  oil,  the  surface  developed  patches  of 
nucleation  sites  that  spread  slowly  with  increasing  heat 
flux,  until  they  covered  the  entire  surface  (point  C1). 
Region  E-E  (or  D’-E)  shows  that  oil  increased  the  wall 
superheat  slightly  for  3  percent  oil  and  significantly  for 
10  percent  oil.  Again,  this  agrees  with  the  concept  of  an 
increased  mass  diffusion  resistance  by  the  addition  of  oil. 

Figure  6.2  shows  the  heat-transfer  coefficient  of  the 
smooth  tube  in  R-1 14-oil  mixtures  as  a  function  of  heat 
flux.  The  curves  show  the  heat-transfer  coefficient  of  the 
smooth  tube  in  region  D-E,  after  complete  nucleate  boiling 
has  been  initiated.  The  effect  of  adding  less  than  6 
percent  oil  is  seen  to  be  small  (about  a  10  percent 


reduction  in  the  heat-transfer  coefficient).  However,  an 
oil  concentration  of  10  percent  causes  a  more  significant 
drop  in  the  heat-transfer  coefficient  (from  0  to  35  percent 
depending  on  the  heat  flux)  . 

Figure  6.3  shows  more  easily  the  degradation  that  oil 
causes  in  the  boiling  heat- transfer  performance  of  the 
smooth  tube.  This  figure  plots  the  heat-transfer  coeffi¬ 
cient  of  the  smooth  tube  relative  to  the  heat-transfer  coef¬ 
ficient  in  pure  R-114  as  a  function  of  oil  concentration. 
The  effect  of  oil  can  be  seen  to  depend  also  on  the  heat 
flux.  Oil  can  be  seen  to  generally  degrade  the  performance 
of  the  smooth  tube,  except  at  heat  fluxes  less  than  5  kfi/m2 
and  oil  concentrations  between  2  and  8  percent.  This 
behavior  was  also  seen  by  Henrici  and  Hesse  (see  Figure 
2.2) .  Since  this  curve  shows  the  heat-transf er  performance 
in  region  D-E,  with  complete  nucleate  boiling,  tne  non¬ 
linear  variation  of  the  physical  properties  of  refrigerant- 
oil  mixtures,  including  surface  tension,  again  probably 
accounts  for  this  anomalous  behavior.  Since  no  measurements 
of  the  physical  properties  of  the  R-1 14-oil  mixture  were 
made  during  this  investigation,  any  possible  non-linear 
property  variations  of  the  mixture  used  are  unknown, 
requiring  future  work. 

B.  BOILIHG  PERFORHAICI  OF  HIGH  FLUX  SURFACE 

Figure  6.4  shows  the  nucleate  pool-boiling  performance 
of  the  High  Flux  surface  in  R-1 14-oil  mixtures.  The  small 
magnitude  of  the  wall  superheats  obtained  during  nucleate 
boiling  should  especially  be  noted.  The  High  Flux  surface 
in  pure  R-114  showed  typical  natural-convection  region  (A-3) 
behavior.  The  incipient  point  (B)  occurred  at  much  lower 
heat  fluxes  and  superheats  than  for  the  smooth  tube  as  shown 
earlier.  The  transition  to  nucleate  boiling  (B-C)  was 


similar  to  the  smooth  tube,  occurring  very  rapidly  (less 
than  1  second)  .  Nucleate  boilng  from  the  High  Flux  surface 
results  in  a  drop  of  the  wall  superheat  by  about  a  factor  of 
10.  This  is  due  to  the  extremely  high  density  of  nucleation 
sites  present  on  the  High  Flux  surface.  From  the  low  heat 
fluxes  (1-10  kW/m2)  at  which  the  transition  to  nucleate 
boiling  occurs,  it  can  be  seen  that  the  High  Flux  coating 
also  assists  in  the  activation  of  stable  nucleation  sites. 

adding  oil  to  the  High  Flux  surface  appears  to  delay  the 
onset  of  nucleate  boiling  (point  B  wall  superheat  increases 
with  increasing  oil  concentration).  However,  the  transition 
to  nucleate  boiling  still  occurred  very  rapidly,  even  at  10 
percent  oil.  The  rapid  transition  to  nucleate  boiling  (less 
than  1  second)  on  the  High  Flux  surface  is  probably  due  to 
the  interconnected  cavities  which  can  assist  in  nucleating 
the  entire  surface  once  a  single  site  becomes  initially 
active.  Oil  is  unlikely  to  inhibit  this  characteristic  of 
the  High  Flux  surface  though  it  apparently  delays  the 
initial  activiation  of  the  first  nucleation  site.  As  seen 
in  region  D-E,  for  heat  fluxes  less  than  37  kH/m2,  the 
effect  of  adding  oil  to  the  High  Flux  surface,  once  it  has 
been  nucleating  fully,  is  to  cause  about  a  30  percent 
increase  in  the  wall  superheat.  At  heat  fluxes  in  excess  of 
37  kw/m2  and  at  10  percent  oil,  the  wall  superheat  increased 
dramatically.  As  seen  in  Figure  6.1,  the  wall  superheat  at 
a  heat  flux  of  98  kff/m2  and  10  percent  oil  is  about  the  same 
for  both  the  smooth  tube  and  High  Flux  surface. 

Figure  6.5  shows  the  heat-transfer  coefficient  of  the 
High  Flux  surface  as  a  function  of  heat  flux  for  various  oil 
concentrations.  Again,  oil  is  seen  to  degrade  the  nucleate 
pool-boiling  heat-transfer  performance. 

Air-conditioning  plants  typically  operate  in  the  heat 
flux  range  of  1 0  to  40  kW/m2  with  less  than  1  percent  oil. 
In  this  region  of  pracical  interest,  the  boiling 


heat-transfer  coefficient  of  the  High  Flux  surface  is  about 
10  tiaes  better  than  a  smooth  tubei  Also,  it  experiences 
only  a  20  percent  drop  in  the  boiling  heat-transfer  coeffi¬ 
cient  with  the  addition  of  oil  (1  to  10  percent).  This  is 
seen  clearly  in  Figure  6.6  which  plots  the  heat-transfer 
coefficient,  relative  to  the  heat-transfer  coefficient  for  0 
percent  oil,  as  a  function  of  oil  concentration.  From  1  to 
10  percent  oil  at  a  heat  flux  of  14  kR/m2,  the  heat-transfer 
coefficient  is  about  80  percent  of  the  no-oil  heat-transfer 
coefficient. 

Figure  6.6  also  shows  that  the  oil-caused  degradation  of 
performance  on  the  High  Flux  surface  is  nearly  independent 
of  oil  concentration  at  practical  heat  fluxes.  Only  at  a 
heat  fluxof  98  kW/m2  and  6-10  percent  oil,  does  the  High 
Flux  surface  experience  significant  performance  degradation. 
At  high  heat  fluxes  and  oil  concentrations,  the  effect  of 
oil  may  be  to  "clog"  the  interconnecting  cavities  of  the 
High  Flux  surface  due  to  boiling  off  of  the  R-114.  Clogging 
the  R-114  surface  with  oil  would  prevent  replenishment  of 
the  nucleation  sites  with  R-114  liquid,  preventing  the 
nucleation  process  and  leading  to  higher  superheats.  The 
time-dependent  behavior  of  the  High  Flux  surface  in  high  oil 
concentrations  and  at  high  heat  fluxes  was  not  studied  in 
this  experiment. 

C.  COMPARISON  OF  HIGH  FLOX  TO  SMOOTH  TO BE  BOILING 

PERFORMANCE 

Figure  6.7  shows  the  heat-transfer  performance  of  both 
the  High  Flux  surface  and  the  smooth  tube  as  a  function  of 
heat  flux.  Again,  the  7-10  times  improvement  in  the  heat- 
transfer  coefficient  by  the  High  Flux  surface  is  easily 
seen.  At  extremely  high  heat  flux  and  high  oil  combinations, 
the  High  Flux  surface  performs  only  sligthly  better  than  the 
smooth  tube. 
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Figure  6.7  also  shows  that  the  heat-transfer  curves  of 
the  High  Flux  surface  and  the  smooth  tube  are  not  parallel. 
At  low  heat  fluxes,  the  High  Flux  surface  is  more  effective, 
in  comparison  to  the  smooth  tube,  in  enhancing  nucleation 
than  at  moderate  heat  fluxes.  At  high  heat  fluxes,  the 
performance  of  the  High  Flux  coating  begins  to  converge 
toward  the  smooth-tube  curve  because  the  surface  became 
vapor  blanketed.  This  is  expected  since  both  tubes  should 
perform  about  the  same  when  they  are  completely  vapor 
blanketed. 

Figure  6.8  shows  the  relative  improvement  of  the  High 
Flux  surface  over  the  smooth  tube  as  a  function  of  oil 
concentration.  For  the  heat  fluxes  of  practical  interest  in 
air-conditioning  plants,  the  High  Flux  surface  is  7-10  times 
better  than  the  smooth  tube. 

0.  EFFECT  OF  SATUBAT20M  TEMPERATURE  OH  BOILING  PERFORMANCE 

As  reported  by  Stephan  [Bef.  9]#  the  effect  of 
increasing  the  saturation  temperature  for  both  the  smooth 
tube  and  High  Flux  surface  was  increased  heat- transfer 
performance.  Figure  6.9  shows  the  improvement  in  heat- 
transfer  performance  in  pure  R-114  achieved  by  raising  the 
saturation  temperature  from  -2.2  °C  (28  <>F)  to  6.7  °C 
(44  OF).  At  high  heat  fluxes,  little  improvement  is  seen  in 
the  High  Flux  surface  performance  because  the  surface  is 
nearly  vapor  blanketed  with  bubbles. 

Figure  6.10  shows  the  effect  of  saturation  temperature 
on  R-114-oil  mixtures.  Increased  saturation  temperature  is 
again  seen  to  improve  the  performance  of  the  High  Flux 
surface  as  well  as  the  smooth  tube,  even  with  10  percent 
oil.  This  is  consistent  with  the  no-oil  results,  but 
contradicts  Henrici  and  Hesse’s  data  (see  Figure  2.5). 


Figure  6.8  Belative  Iaproveaeut  of  Higl 
Surface  over  Saooth  Tube  versus  Oil. 


E.  10CAI  OIL  SAHPLI1G  EFFORTS 


The  results  of  the  attempt  to  sample  oil  locally  near 
the  boiling  tube  were  inconclusive.  The  method  outlined  in 
Section  IV. C  did  not  produce  either  repeatable  or  accurate 
results.  Checks  of  the  sampling  method  were  made  by 
sampling  the  bulk  liguid  with  0  and  10  percent  oil  at  zero 
heat  flux.  The  0  percent  oil  check  yielded  oil  concentra¬ 
tions  from  0  to  2  percent.  The  10  percent  oil  check  yielded 
oil  sample  concentrations  from  5  to  25  percent. 

Further  refinements  to  the  sampling  apparatus  are  needed 
to  permit  an  accurate  local  sample  of  oil  in  the  vicinity  of 
a  boiling  tube.  Section  VII.  B  contains  reconrmendatiors  to 
improve  the  sampling  apparatus. 
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VII.  C08CI05IQNS  AND  RECOMMENDATIONS 


CCICIUSIOMS 

1.  In  pure  R-114,  the  pool-boiling  heat-transfer  coeffi¬ 
cient  of  the  High  Flux  surface  is  about  10  times 
larger  than  that  of  a  smooth  tube. 

2.  The  High  Flux  surface  began  nucleate  boiling  at  low 
heat  fluxes,  about  1  kW/m2,  compared  to  a  heat  flux 
of  about  10  kR/m2  for  the  smooth  tube. 

3.  Oil  delayed  the  onset  of  nucleate  boiling  with  the 
High  Flux  surface. 

4.  Oil  resulted  in  about  a  20  percent  reduction  in  the 
heat-transfer  coefficient  of  the  High  Flux  surface 
for  most  practical  heat  flux  (less  than  37  kw/m2)  and 
oil  combinations  (less  than  6  percent)  . 

5.  At  heat  fluxes  of  98  kW/m2  and  greater  than  6  percent 
oil,  the  performance  of  the  High  Flux  surface 
degraded  by  as  much  as  80  percent.  The  performance 
of  the  High  Flux  surface  was  little  better  than  the 
smooth  tube  at  98  kN/m2  and  10  percent  oil. 

6.  The  boiling  heat-transfer  coefficient  of  the  High 

Flux  surface  was  about  7  times  better  than  that  of  a 

smooth  tube  for  oil  concentrations  from  1  to  10 

percent  over  the  range  of  heat  fluxes  employed  in 
air-conditioning  plants  (10-40  kw/m2) . 

7.  The  boiling  heat-transfer  coefficient  of  both  the 

High  Flux  surface  and  the  smooth  tube  increased  with 
increasing  saturation  temperature.  The  improvement 
decreased  at  high  heat  fluxes  for  the  High  Flux 
surface  due  to  vapor  blanketing  of  the  surface. 


.  Local  oil  sample  results  were  not  accurate  or  repea¬ 
table  enough  to  determine  if  a  variation  in  the  local 
oil  concentration  occurs  near  a  boiling  tube. 

BECCHHEIDATIOIS 

.  The  present  cartridge  heaters  should  be  replaced  with 
more  reliable  heaters  that  will  produce  a  uniform 
heat  flux  axially  at  all  heat  loads. 

.  A  solid  copper  tube  should  be  coated  with  High  Flux 
and  tested.  Eue  to  the  very  small  superheats  of  the 
High  Flux  surface,  even  a  small  amount  of  contact 
resistance  could  affect  the  data  at  low  heat  fluxes. 

.  The  physical  properties  of  the  R-1 14-oil  mixtures 
tested  should  be  measured  to  obtain  information  to 
better  explain  the  reasons  for  the  heat-transfer 
performance  of  both  the  High  Flux  surface  and  the 
smooth  tube.  Particularly,  the  anomalous  rise  in 
heat-transfer  performance  of  smooth  tubes  when  1  to  6 
percent  oil  is  present  in  R-114  should  be  studied. 

.  A  secondary  heater  should  be  installed  in  the  boiling 
section  to  keep  the  total  heat  input  constant.  The 
heat  input  to  the  secondary  heater  should  be 
increased  when  the  heat  input  to  the  boiling  tube  is 
decreased  and  vice  versa.  This  modification  will 
maintain  a  constant  heat  load  on  the  condenser  and 
eliminate  the  need  to  operate  valve  VC,  except  to  set 
the  saturation  temperature  at  the  beginning  of  a  run. 
Alternately,  the  effort  to  computer  control  valve  VC 
should  continue  by  obtaining  an  accurate  thermistor 
with  a  high  temperature  resolution  for  use  with  the 
STS- 2 A  microcomputer.  Improved  temperature  resolu¬ 

tion  would  allow  the  computer-controlled  valve  to 
operate  properly  and  free  the  operator  from  the 
demands  of  manual  control  of  valve  VC. 


The  performance  of  a  bundle  of  High  Plux- coated  tubes 
in  E-1 14-oil  mixtures  should  be  studied. 

Data  should  be  obtained  on  the  High  Flux  surface  over 
a  wider  range  of  temperatures  and  with  oils  of 
varying  viscosities. 

The  High  Flux  surface  should  be  tested  for  time- 
dependent  heat-transfer  performance  by  taking  data 
periodically  over  a  long  time  during  boiling. 

A  small  oil  sample  container  with  valves  and  a  vent 
line  should  be  manufactured.  The  flexible  tubing  and 
pinch  clamps  used  in  this  experiment  did  not  properly 
hold  or  vent  the  sample. 


APPENDIX  A 

THEBNOCOOPLE  CALIBRATION 


Karasabun  £Bef.  8]  describes  the  thermocouple  calibra¬ 
tion  equipment  in  detail.  Two  thermocouples  were  calibrated. 
One  was  made  from  wire  at  the  beginning  of  the  roll,  the 
other  from  the  end  of  the  roll,  following  the  making  of  all 
thermocouples  used  in  the  apparatus  and  tests. 

Essentially,  the  manufacturer-supplied  calibration  equa¬ 
tion  for  the  thermocouple  wire,  a  seventh  order  polynomial, 
was  corrected  slightly  by  adding  to  it  a  second  order  curve 
fit  of  the  variation  of  the  manufacturer’s  computed  tempera¬ 
ture  for  a  given  emf  from  a  known  set  of  reference  tempera¬ 
tures  (measured  using  a  Hewlett-Packard  2804A  quartz 
thermometer  with  a  temperature  resolution  of  ±  0.0001  K  and 
accuracy  of  ±  0.03  K) . 

The  manufacturer’s  emf  to  temperature  conversion  equa¬ 
tion  is : 
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(A.1) 


Figure  A.1  shows  the  guartz  thermometer  reading  minus 
the  thermocouple  readings  (discrepancy)  versus  temperature. 
The  two  thermocouples  agreed  to  within  0.05  K  of  each  other 
and  the  manufacturer's  seventh-order  polynomial  needed  about 
a  0.  1  K  increase  to  more  accurately  convert  the  emf's  to  the 
true  temperature.  The  correcting  second-order  polynomial 


DCP  =  bQ  ♦  b1T  ♦  b2T*  (A. 2) 

where 

DCP  =  discrepancy  (K) 

=  8.6268968E-2 
b  =  3.7619902E-3 

b  =  -5. 0689259E-5 

2 

T  =  thermocouple  reading 

(from  eguation  A.  1)  (<>C) 

Thus,  the  temperatures  computed  by  the  data-reduction 
program  (DRP2)  were  emf's  converted  to  temperature  by  egua¬ 
tion  A.1  with  corrections  for  that  temperature  computed  by 
eguation  A. 2  added  to  the  temperature  to  get  the  true 
temperature. 

Since  the  data-reduction  program  utilized  differences 
between  thermocouples  in  all  computations,  such  as  wall 
temperature  minus  saturation  temperature,  the  corrections 
above  were  necessary  only  for  the  computation  of  items 
dependent  on  the  absolute  temperature,  like  the  fluid  prop¬ 
erties.  Appendix  E  describes  in  detail  the  uncertainty 
analysis  and  the  effect  of  wall  temperature  variation  on  the 
computation  of  the  heat-transfer  coefficient.  Thermocouple 


Figure  A.l  Thermocouple  Discrepancy  Correct 


variations  of  the  soft-soldered  tubes  were  most  likely  due 
to  a  slight  amount  of  contact  resistance  and  to  the  surface 
characteristics  of  the  tubes  tested,  since  data  runs 
involving  shifting  the  wall  thermocouples  did  not  affect  the 
data. 


APPENDIX  B 

DATA  REDUCTION  PROG BAB 

The  data  reduction  program  below  consists  of  the 
following  sections: 

Bain  Program  -  Menu  of  subprogram  options 

Sub  Main  -  Take  data  or  reprocess  data 

Sub  Plot  -  Plot  data  on  log-log  scale 

Sub  Poly  -  Compute  least-sguares  curve  fit  of  data 

Sub  Flin  -  Plot  data  on  linear-linear  scale 

Sub  Stats  -  Compute  average  and  standard  deviation  of  data 

Subprogram  Bain  consists  of  the  following  steps: 

1.  Create  data  file  for  data  and  plot  points. 

2.  Select  tube  type. 

3.  Monitor  heat  flux  or  saturation  temperature  to  estab¬ 
lish  steady-state  conditions. 

4.  Scan  all  channels  listed  in  Table  2  and  save  in  data 
file. 

5.  Convert  raw  emf's  to  temperatures,  current,  and 
voltage. 

6.  Compute  the  heat-transfer  rate  for  the  cartridge 
hater. 

7.  Compute  the  average  wall  temperature,  the  wall  super¬ 
heat,  and  the  film  temperature. 

8.  Compute  the  physical  properties  of  R-114  using  given 
correlations  at  film  temperature. 

9.  Compute  the  natural-convection  heat-transfer  coeffi¬ 
cient  of  R-114  for  the  non- boiling  ends  of  the  test 
tube. 

10.  Compute  the  heat  loss  from  the  non-boiling  ends. 


11.  Calculate  the  corrected  heat  flux  from  the  test  tube 
to  the  liquid  B-114. 

12.  Calculate  the  boiling  heat-transfer  coefficient  of 
the  B-114  fron  the  test  tube. 

13.  Print  data.  Store  heat  flux  and  vail  superheat  values 
in  plot  file. 

The  following  is  a  listing  of  the  complete  data  reduc¬ 
tion  program  (DBP2)  written  in  Basic  3.0  for  the 
Hewlett-Packard  9826  computer. 
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1995 

«_  ■ 

:ooo 

next  Ji 

:oo5 

IF  I-l  THEN  Volt»5ii««5 

V. 

:0«0 

If  I»?  THf N  Anp^C 

»  * 

2015 

NEXT  1 

■-  ’ 

:02  c 

Aodo-Volt  •Ano/T  r  I  •l)2,L  1 

1 

2025 

TF  A8S< Aadn-Oado >  Err  THEN 

2030 

If  Aqdo  Oodp  THEN 

’  ' 

2035 

BEEP  4000  .2 

.. 

2040 

BEEP  40P0  -2 

2045 

BEEP  4000  .2 

:os0 

ELSE 

. 

2055 

BEEP  250  2 

2060 

PEEP  250,  .2 

2055 

BEEP  2S0..2 

r 

2070 

ENP  IF 

A 

2075 

PRIN!  IJSIMG  *4X  H2.30F.2M 

. 

2090 

WAIT  2 

2095 

GOTO  1970 

2090 

ELSE 

2095 

PEEP 

2100 

PRINT  USING  "lX.HJ.3Pf. -C*  . 

2105 

Err»SP0 

2H0 

WAIT  2 

21  »5 

GOTO  1970 

A 

2120 

CNO  IF 

*  • 

2125 

2130' 

END  IF 

w. 

2  f  -•  5  • 

LOCI'  10  c€  J  Ti*t 

2140 

i r  ido-:  then 

2U5 

BEEP 

2150 

INPUT  TNJtR  OESIRED  Tsai 

2 1 55 

PRINT  USING  DT»at 

• 

s 

r  ico 

Oldl '0 

2155 

0 1 1?2  *  0 

2170 

OUTPUT  709 i ‘ AM  AE33  M  35  1 

2 1  75 

FOR  1*1  TO  3 

V- 

2 1  90 

5-1W--0 

•  V 
fc1 

2ie5 

OUTPUT  709 « 'AS  SA¬ 

2130 

POR  I i-1  TO  2 

2155 

ENTER  7P9 1 E  lio 

2200 

Sort"5u«*E l to 

2205 

NEXT  j, 

U  2:10 

E  l  ta-fjnn/2 

p 

22)5 

T 1  d*ENl  vjv<  F  1  i  o  » 

r 

2220 

If  |«1  IHFfl  Af  Id-1  Id 

r- 

2225 

IE  I *3  THf N  Tv ■ T 1 d 

r-  :::o 

NEXT  I 

i*. 

22  35 

IF  AP«<  At  Id-fU  Id  •  •  .2  THf  N 

•T. 

2240 

IE  AIM  Ot  Id  THEN 

2245 

Pf f  E  4000  .  .1 

2250 

BEEP  4000 

— 

2255 

TFEP  4000 . . 2 

P 

3250 

Lf  T 

2265 

ec'tp  250  .2 

GFEP  250,.: 

’  “  ,'5 

pi  rr  :rr>  ,r 

2  ~r0 

Ff«0  If 

•* 

-  -  {'i; 

E  *  r 1 • A  t  1 d  Oldl 

i 


Pate 


'US' 


v\ 
"  -y 


e f r : * f v  ni.12 

::30  oi-im* 

pnrnr  umng  'J»  s-hoodo  o»">ond  miu  e»ti  .n  Err.. 
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WAIT  : 

G.ifO  2*  70 

else 

tF  ABS< At  Id  Ot Id  I  .1  THFN 

IF  At !d-0* id  THEM 

BEEP  3000..  2 

BEEP  3000  .2 

ELSE 

BFE.P  000,  .2 

trcp  soft,.: 

END  IF 

Err  I  “At  Id  -0  Leif 
Oldl-At  id 
Err:*Tv-0ld: 

OldZ'W 

PRINT  USING  "  4X , 5<  MOD . 00 . 1  *  ) “ i Dk 1 d , A 1 1 d  .Err  I  .  Tv  .Err 2 

uaii  : 

GOTO  :i70 

FISE 

BEEP 

Err l -At id-Oldl 
Ol.ll-Mld 

Fr-:-Tv-oiu:  * 

c  i  r  v 

PRINT  USING  *4X  .5'.  Mnn.f.'O  tX)' iOt  H  AMd  Errl.Tv.F-r2 

VAIT  2 

0010  2170 

ENP  IF 

ENO  IF 

END  IF 

FPFOP  TRAP  FOR  I.jo  OUt  Of  PO"NOS 

rr  i .jo  :  thin 

PEEP 

r.r.tn  ipq® 

F no  ir 


2  4^01 

I#iKF  PA  1 A  ir  In 

-0  l 

2J9S 

pfet 

27-00 

INPUT 

FoirR  r>»lk  o 

:s05 

m.urin 

70S i 'AR 

af:s 

:sio 

f  OR  l  - 

1  If)  tZ 

:?is 

OUTru? 

70S i  AS 

SA- 

:s:o 

Son ■=  0 

:s:s 

ron  Ji 

-i  to  : 

2ST0 

f.NTFR 

70S.E 

:s:s 

Sun  *F.infE 

:S40 

MF<t  ? 

1 

Fnf  (  I 

•Sum/: 

nr  <  r  i 

27^S 

OUTPUT 

70s,  sn 

APS: 

:sfc 

FOR  1  * 

i  rp  : 

rsr.s 

OUTPUT 

700 s  AS 

S  A“ 

:s-0 

*  lJf»‘0 

•'S'S 

rt’R  le 

- 1  to: 

2SF0 

ENTER 

”?9iE 

7  7  PS 

Son-  on*F 

2r  *'0 

•if  v  I  .1 

» 

?Sir, 

[P  1*1 

r\H  i*r- 

r:-m 

:c00 

ir  !>; 

TUFN  Ir- 

A.jf, 

'HOS 

ur  a  r  i 

l  Bon 
AL  3R  UR S' 


ClIfCR  AT  !  J  -  ;  Rr,n  T  ■  -  |  Fl*tfi 

run  ir 


CT"Tni  «(irc  tu  ir  mi-'  ■•pi  i  .'I'hprur 
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:640 

FOR  I  *  1  TO  17 

:64s 

IF  I d* c  0  THEN 

:6so 

IF  I  -|.dtc  1  OR  I-l.dlc: 

2655 

T (  I  »*- 59. 99 

2660 

6010  2710 

2665 

END  IF 

2670 

END  IF 

2675 

IF  lit  4  THEN 

26B0 

IF  1)4  AMO  19  THEN 

2$R5 

T (  I  >*-99.99 

2610 

GOTO  2710 

26  9S 

EMO  IF 

2700 

END  IF 

2705 

T<  I  )-FNTviv<  Fnf ( I  )  > 

2710 

NEXT  I 

2"  1 5 

IF  Itt<4  THEN 

2720 

FOR  1*1  TO  4 

2725 

IF  l-Ldtcl  OR  l-ldtc 

27^0 

T  hi£*Tw<) 

2775 

else 

2740 

T  wa* Tud*  T i  I ’ 

2745 

END  Jfc 

27C0 

NFV  T  I 

:rss 

T  na  *  I  wd  /  (  4-Idtc  ' 

2760 

n  ;r 

2765 

FOR  1*1  ro  5 

2770 

IF  I  «t.dtp  1  OR  l-Ldtc 

2775 

T  ud  *  T  b>d 

2-40 

ri  5C 

2775 

r..a-=fw-*»f<  i  * 

2-90 

F NO  IF 

2715 

nr  <i  i 

7?0P 

TM-'TrtMfl  Ixlte  » 

2  ?  25 

END  IF 

^010 

TH*rn  > 

2  P  1 5 

T  M2  T :  1 0  * 

;?n_*0 

f^TMl  • 

2025 

IF  Itt  l  IHIN 

25  20 

Tld7*- 99.99 

2?  ?5 

T  »(  It  I0'*T<  II  i'7 

7<?40 

END  IF 

7745 

T  *»tmn  TI  1  2  * 

2-50 

Amo  *  I » 

:?i5 

1  t  -'»r  .  25 

"550 

■  l  t  «  Amo 

2P65 

IT  1 1  t • 0  THFN 

2970 

V r  ij*f  NK  C'.l '  1  l»M  > 

:o,5 

Ft  SF 

7590 

F'  r; u  **  “*»dt  !  *  l  ' 

2205 

F NO  IF 

:^?i 

2 ‘*'36  i 

mttpiFR  i’^MD'.'CMOll  1 

2 

U.«  1. .a  £••  Of-  07  01  - 

2T75 

I  hct  v'h  •  r  ..  1  1  1 

Input  MAIL  I A I  •  ifl r*»*  T  J  ril  »p  I  -TNI?  - 

lr  rum  to  !*• 

If  t-.  i  u  if  w  ’“*0 
tun  ip 

r'*np*t|E  ''An 1 011*5  fcnriRtic? 

»  f  I  |m  -f  nr  f  i  I  ml  1  M  f  hi  ' 

Rli?  ‘  Ml- hr  ••  f  f  I  I  *  1 
M»i-  fhHm  T  f  i  I  n  1 
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THE  INFLUENCE  OF  OIL  CONTAMINATION  ON  THE  NUCLEATE 
POOL-BOILING  BEHAVIOR  0.  .  (U)  NAVAL  POSTGRADUATE  SCHOOL 
MONTEREV  CA  J  T  REILLV  MAR  85  NPS69-85-002 

2/2 

UNCLASSIFIED 

F/G  13/1 

NL 

2970  Co-FNCoi Tf »1«> 

2375  Bete-FNeetalTfiJw) 

29E0  Ni-Mu/Rho 
299S  Aloha-K '<Rho»Co > 

2<190  Fr-Ni/AU-hd 
2995  P*a**FNPaat<  T!d> 

5000 « 

3005'  COMPUTE  NATURAL -CONVECTH'G  HEAT- rRANSFER  COEFFICIENT 
3010'  FOR  UNENHANCED  ENO<SI 
3015  Mbar-19# 

3020  Fe*<  Hbar*P*  *  Kcu* A  )  ?  * .  5«t.u 
3025  T anh-FNT anh( Fc  > 

3030  IhetaMbelab^Tarm'Fe 

3035  V  •  •<  9.81  *0el  *?•  thetab'Oo  5*  f  arth/ i  Ee*Ni  'Aloha  .• »  .I6R667 

3040  Vy»<  I  ♦€  .559/ Pr  >  i  9MB  » )  (  8  27  I 

3045  Hb*rc*K'0e-«  .R*.387«X.-  'rv  )  2 

3050  IF  APS(  <  Hba-'  Htoarc  )/Hbarc  )'  .001  THEN 

3055  Hbar«CHbdr«Hharc »* .5 

3050  GOTO  3020 

3055  ENO  IF 

3070' 

3075'  COMPUTE  HF.*T  1.055  RATE  THROUGH  ONE  NUANCE  D  END'S? 

2PC0  QI-<Hhar*P«Kcu*A >  .S«Thetab« lanh 
3005  Oc-Q  -2*01 
3090  Aa*°I *02  *L 
3P9S* 

3100'  COMPUTE  ACTUAL  MEA1  ElU*  AND  BOILING  COEFFICIENT 
3105  Qdn-Oc  -'As 
3110  Htuhe-Odo/ Thetab 
31  15' 

n:0»  pecoro  riME  or  data  iamng 

3125  IF  Irt-P  T  NFN 
3130  OUIFUT  .0^1 "TO* 

3 1 "5  EMIER  709* I o 13* 

31 40  FUO  n 
3U5' 

3150'  ourrur  oata  r«j  prinifr 
3lr-5  FRINIER  IS  701 

3i c-e  if  i-v  0  mrri 

.-165  rut  til 

x 1 70  PRINT  vising  M0«.*“Oa*a  Set  tlunber  -  “ '  ,000  2X.**Bull  Oil  X  ■  "  “  0O.D.5X.I 
in'  » I  Bob 
5*75  PRINT 

3 1  P0  PNINf  USING  *  I  0 K  ,  *  “  T C  flo  :  1  2  3  4  5  5  7 

9  “  “  ' 

3195  PRINT  USING  "  1  0*  .  “  *  I  cmu  8C  t  x  .MPD .  00  ’•*  i  T<  I  )  .  T(  2  >  .T<  3  > .  T  <  4  )  ,  T(  5 '  .  T<  6  >  ,T  * 
7  •  ,T  .  P  ) 

’IIO  PRINT  U*.:IN»5  'fOV  “ "  U».-  Uiod  1  !  iud2  I*  jo r  FJsat  T^unp**' 

1R5  PRINT  USING  "  IPX  Mf’O  .I'P  U  l  ,IX  .MT‘0  .  GO  .  1  X  .2  *  IX  ,MP0 . DO  '  ,Zx  .  MOD  .0"*Twa.T|d, 

I  IdZ  ,  Tv  ,P'.  ji  T  viao 

3200  PRINT  USING  “  1 0X  •-  Ihi-t^b  Mtube  Qdo"  ' 

3205  KP  HIT  USING  10*  1100.70  IX  .  M  Z .  3PE  .  I '  . H J  .  30E  '  i  T he ♦  ah  Htube  Odo 

32  1  0  l  NO  11 

3215  IF  Jc-"l  IhFN 
3230  IF  .1-1  THEN 
32.35  PfJlfll 

3230  PPINI  USING  ”I0X  '"RUM  No  OilX  Tsat  Htube  Qdo  Thelab*" 

33  35  EMP  II 

7340  FPINT  USING  •I7.X  !0  4*  DL1.2X  MDO.IiD  .  3t»X  ,MZ  .  30E i  J  Boo  Md  Htube  Qdo. That 
ab 

3315  CNO  U 

r:r«  if  j«  0  then 
::rr.  pf  ep 

32FO  I hF" I  ' OK  fn  5  ! ORE  r n f r  DATA  *E  f  <t--f  0 -N  <->» 

3255  END  If 

7:  /0  IF  0 1  -  1  OH  l«i  I  ff'T  N  ! 4  1  H 
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3275  IF  PT»l  ANO  ln*0  THEN  (Hi  1  PUT  9FiletiGoo  Tcldt  E*f(  •  t  ;Vr  ,  Ir 

3390  IF  In- 1  PR  0*  *  I  THEN  OUIFIJT  •Plot  iQdn  :  Thet  ab 
3205  IF  ln*0  THEM 
3390  BEEP 

3395  I  Fir  U 1  "WILL  THERE  BE  ANOTHER  RUN  <  1  -  Y  .  0«N  )’»  “  .Go._on 
32C0  Nrt»n*J 

3305  IF  Go. on*  I  THEN  3330 
3310  IF  6o_on-1  THEN  Reoeat 
3215  ELSE 

3330  IF  J.NrunM  I  HEN  Rente t 
3325  ENO  IF 
3230  IF  l«*0  THEN 
3335  BCEP 

3340  PPinr  USING  "1 04. ""NOTE:  "-.27,*"  data  funs  were  stored  in  rile 
I ,02. f i lef 

ASSIGN  «rilel  TO  • 

0U1 FUT  #f i le2»Mrun-l 
ASSIGN  9f  1 1  e I  TO  01.  I  •. 

ENTER  fF»  lei  tO.it-*  .1  ltd  Idtcr.Itl 
OUTPUT  Af»1rr»Pa»p*  I  dtr  I  ,1  dl«2 . 1 1 1 
FOR  1*1  TO  Nr un- 1 

ENTER  <*Filtfl.‘Bno  TolU*  ?«/<•>  Ur  Jr 
pti!f»n»  0fi|r2tPoO  Tolrlf  rmf  <•>  ,l,r  Jr 
N  fVT  I 

ASSIGN  0T  i  lei  TO  • 
pupr.£  -DUMMY ' 

ENO  IF 

Ltd  ,  4 

r-R  i  m  r  ijSING  '  1 0X  ,  *  “NO  IE  J  "  *  .11  ,  were  stored  in  plot 

••  "  10A"  i  J~  \  P_f\  le* 

ii  is  A-;e.|r,N  <»r  ,  je?  in  • 

3J20  ASSIGN  pilot  10  « 

3425  CAl  l 
3430  PFEP 

34:5  INPUT  I  1*F  TO  PI  OT  DATA  (1  *Y  ,0-N  » > "  ,01 
3440  IT  01  *1  IMFN 
3445  CALL  Ptct 

J4'fl  EMp  |F 
3455  ci|HfNl) 

3460  ■ 

?4F5i  CUH"t  MIS  Of  PROTTRlV  P'NlJIoNS 
2470  off  rwcut  i ' 

24 ’5'  OF HC  COPPER  250  TO  200  ► 

34E0  Tl  1*273. 15  T  TO  ► 

3405  h  *4/4  .112*11 
34R0  RETURN  K 
2495  FNf NO 
3500  OFF  T  Nf1u<  T  > 

?505»  170  10  ;R0  k  curve  in  ui  ui’irnurjN 

3510  TIM *273. 15  'C  TP  K 

3515  Ng»E  *P<  -4 . 4G36*<  IP*  I  .  47  ■  Tl  *#»•». Of  - 3 

3520  PFll'RN  Mu 

7525  I  Nr HP 

3530  off  rnro< 1  * 

35351  11*0  TO  400  V  ritP-'E  Ml  OF  Co 
35  40  TIM  1 273.  «*♦  1  f  I  o  F 

554c;  Co-.  40IJ?«»  t  -B500TF  3*  M  '  1  .514945  -5*  Tl  2 -fi  .  67^76  - 1  0»  T I  3 

7r’r-0  Lp  - 1. 1* '  I  000 

3555  PI  TURN  >  0 

7550  I  NF  HP 

-5fiS  Off  ft  iRI.nf  f  > 

rr  0  ri  1  •:?.?. i«5  ,r  to  ► 

••*375  X-I  •  I  .  0 •  f I  753.95*  *F  »•?  R 

35«'0  Rr  T,.  32*61 .14MI4**  <»  7  t  UR  .  4 1  oOlf.x  *  1 7 . 4. -Pr7R*«  5H.MMP20 
•r^jc;  p,  Ur,  .  .  (3534  ;R 

75^0  Rf  TURN  Ro 


3J45 

37-50 

3355 

3260 

3365 

3370 

3375 

3700 

33F5 

37.nO 

3795 

3400 

3405 

2410 
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3595 

3600 

3605 

36 1C 

36 1 5 

3630 

3625 * 

3G30 

3675 

7640 

364*5 

3650 

3665 

3660 

3665 

3670 

3675 

3EP0 

36*5 

3590 

3635 

3700 

3705 

3710 

37  ?  5 

7770 

3725 

3730 

37^5 

3740 

3745 

3750 

3755 

3760 

3765 

377C 

3775 

7"90 

37P5 

3790 

37951 

7600 

3905 

3010 

39 1 5 

3620 

3925 

3P70 

2975 

3040 

3645 

3?r»0 

3965 

7PG0 

’,"65 

7P70 

<975 

3690 

39  90 
<995 
-.500 

<905 

3310 

■'Siq 

<920 


FNENO 

DEF  FflPrl  f  » 

Pr *F NCp<  I  >*FNMu<  T  >/FNK<  T  J 
return  Pr 
6  ME  NO 
DEF  FNKtT) 

T  360  K  WfJH  T  IN  C 

K- .071-. 000261 •! 

RETURN  k 
FNE  NO 

OEF  FNTanh\X > 

P-EX P< X ) 

Q-l/P 

Ttfnh-«  P-Q)/<P*0  ) 

RETURN  Tanh 
FN6N0 

OFF  ENT/s.<V> 

COM  /C c/  C<7  > ,  (cal 
r-t  ( C  > 

FOR  J«l  TO  7 
T-T*C< I >«W ' l 
NEXT  I 

!F  !c*l»l  THEN 

!-r-6.74?:534E-2fT*(9.0777043E-3  -T*i  9.32539I7E  5) > 

ELSE 

I »T . 626897E-2+T •<  3 .76 1 99E- 3-T*5. 06R9259E -5 > 

ENO  IF 
RETURN  f 
FMfMO 

OEF  FNBettfi r  5 
ftcP«FNftho( T* . 1 > 

Rom'»fhri»o(  t-.  1  > 

Beta-  2/< Rco»Ru* >•< Pcu -Ron >/ . 2 

RFT"RN  Beta 
F  ME  NO 

OFF  Fiirft  Ini  Fu.lld) 

Tf  1  l  ««<  TwlMtft '? 

RETURN  rrti« 

FHFM» 

OFF  FNPiatt  Tc  > 

0  TO  00  deg  f  CURVE  FIT  Of  P**t 
rr-l J2 

Pa-S.  945525 ♦  T f  •<  .  153520P3M  f  •<  1  .  4B40Mr>.5E-3t  ?  f  *9  615067 IE-6  )  ) 
Pg-P*-14 .7 

IF  Pt 1  0  THEN  I  «>P6lG,-*irt  Hq 

PsaT-Pg 

FL  rF. 

P**» -Pg«2 9.92  14.7 
END  If 
RETURN  °sat 
EMEND 
51 W  °lnt 

Pf'M  'Cn!*'  A<10  101  C<  10i  fil  4  >  ,Nop  .  lorn  |  .  Ogo  Hog 
INTfr-FR  Jj 
PRINTER  15  1 
en  p 

IfJPlir  "I  I K  F  OFF  AMI  T  UAlilf'.  F‘>P  PI  01  (  »  •  f  ,  0-N  >7-  \  rfv 

eEm 

PRINT  Ilf, IMR  -4x  “C-l-rf  n,|||i*ni"' 

TPINI  l*c  1  HP  6X  **0  (i  vervn  delta-T""" 

PR  INI  1  '5  INC?  6x  '‘1  h  ^e»  ^ts  deita-T*"* 

n?im  UFTNC.  " 6 V  *  * 2  h  .*».„*  n“' 

ItJpUT  Opn 

PFl  P 

INF"  I  ‘'RFC!  i  IN  I T  6  1  0-5 1  .  >  «H4fil  !  *M  i-  (on 
PRINTER  1?  7 05 
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393*  BEEP 

3935  INPUT  'ENTER  NUMBER  OF  CYCl  ES  FOR  x-A*I5" 
394C  BFEP 

3945  INPUT  “ENTFR  NUMBER  OF  C/OlfS  FOR  Y-AXIS- 
3950  PEEP 

39F5  rUPUT  'ENTER  MIN  X-i.ifll.UF  <  MULTIPLE  OF  IQ) 
3%0  PFFP 

3965  INPUT  “ENTER  MIN  Y-UALUF  <MU1TIPIE  OF  IQ)‘ 

3970  EL9E 

39/5  IF  Qoo-0  THEN 

39P0  Cy-3 

3995  C*-3 

3990  X«»n-.| 

3935  Yi*»n«10® 

4000  ENP  TF 
4005  IF  Ono-1  THEN 
40 1 0  Cy-3 
401 5  Ca-3 
4030  X«in» . I 
4035  Y«i  »»■  1 00 
40  '0  ENP  IF 
4035  IF  Opo»:  THEN 
4040  Cv»3 

4045  c*»r 

4050  Y m i m *  I  000 

4  055  YiMn-100 
40E0  ENP  IF 
1065  EftO  IF 
4070  PFFP 

4075  PRINT  "INiSPliIP  '’300  7700  .8300 .59001  * 

40P0  PRINT  'SC  0,100  0  ,  I00i  Tl  ?  .0|* 

4005  5f  •  » I 00 'C  * 

4090  cfy'100/Cy 
4095  PRINT  "PU  0.0  P0“ 

4100  Nn-9 

4105  FOR  | .1  in  r. M 
4110  X4»*xni„,f0f  I- |  | 

4115  IF  T-r .  H  ri»FN  Nn-» 

4  1 70  FOR  .1-1  TP  fin 
4 1 75  IF  |-|  MIEN  PRINT  “  ff  .7  0" 

4130  JF  .»-?  THEN  PRINT  'TI  1  0- 

4  135  f.l 

4140  X»|  OH  Ka/Vnin  >*SE  < 

4145  PRINT  'Pfl' jX 0i  X T  i  * 

4160  nr V  T  J 
4155  NEXT  I 

4160  PRINT  ‘Pfl  100  fliPtli" 

4165  PRINT  'PU  Pfl  0  0  Pp- 
4170  Nn-9 

4175  FOR  I-i  TO  C»  M 
4IP0  Y<-  *  -  Yr*  i  n»  I  0  '  f  T  -  I  ) 

4195  TF  I  -P v  ►  I  THEN  Nn-| 

4190  FOP  |-f  TO  Nn 
4  195  IF  l-l  TllFN  PRINT  -  Tl  7  0' 

4700  IT  17  IHf  N  ERJMT  -||  I  fl- 

4705  • f 

4 ?  I  0  v-t  RT(  Ya/Vmln  * •  5  f  y 
4715  F R ! M T  *Pfl  0  ' iV  'YT* 

4770  NEXT  .f 

4775  NEXT  l 

4V0  PRINT  “PA  0  100  ri  0  .7' 

4  7  75  Nn  *  9 

4740  FOP  1-f  TO  T. • I 
4745  *M'Mi»in*l0  (  I  -  1  » 

4760  IF  I-P.4I  THEN  Nn-I 


,G.< 

Cy 

’  ,Xntn 
1  .Yntn 
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4355  FOP  J-1  TO  Nn 

47fi0  IF  J-l  THFN  PRINT  'Tl  ?  r 

4755  IK  S'  I  THEN  PRINT  *TL  0  1“ 

4270  Xa«Xat*.l 

4275  X«l  GTf  x«/y*ln>«Sf» 

4200  PRINT  “PR-iX,* ,1001  XT* 

4205  NF*T  I 
4290  NEXT  I 

4295  PRINT  *PA  100.100  P»J  PA  100  0  PD' 

4300  Nn«9 

4305  FOR  1-1  TO  Cy»l 

4.310  Ya*»y«ln*l0*<r-M 

43*5  IF  t-Cy*1  THEN  Nn-I 

4320  FOR  J-l  TO  Nn 

4375  IF  J-l  THFN  PRINT  * TL  0  2" 

4230  IF  J>1  THFN  PRINT  "TL  0  1" 

4335  Y.s*vdt«J 

4340  Y-l  GTf  Ya/Ymin  »*Sfy 

4345  PfliNT  'PO  FA  100,"  i.d" 

4350  NFXT  J 
4255  NFXT  I 

4350  PRINT  'PA  100.100  PIT 
4  Mi 5  rPIMI  "PA  0  *2  5R  I  .5  7" 

4370  J  i  *1  GT<  X«,  n  ) 

4375  FnP  l-l  TO  C- ♦ I 
47P0  Ya'Xnin* 10  '< II) 

4395  X*tGT<Xa/Xmn)«Sf> 

4390  PRINT  *PA' iX  0t " 

4335  IF  Ii'-0  THEN  PRINT  TP  -  2  ,-2  ;L.Bl0iPR  -2  ,2  il.R"  1 1 1  j  “  * 

4400  IF  Ii.0  THFN  PRINT  TP  -  2  . -2  il  R 1 0 1  PR  0  ,2  il  P"  1 1 1 1  * " 

4405  11-liM 

4410  NEXT  1 

4415  PRINT  "PH  PA  0.0" 

4470  I i »l PT< Ynin) 

44?5  Y10-10 

44^0  FOR  I-I  TO  Py  +  1 

4415  Y*-Yi«tn«  10  I  I  -1  > 

4440  VtRT<  Ya/Y«in  )*5fv 
4445  PRINT  'PA  0 , " l V , “ " 

4490  PRINT  "CP  -4  ,-.2Sil  R1CUFR  -2  ,2 iLP* » 1 1 » " " 

4455  l i  - 1 1 ♦ I 
44K0  NEXT  I 
44155  RFFP 

4470  INPUT  "UANI  U3E  OFFAHI  T  I  ABFl  5  (I-Y  0»N>9".fdl 
4475  IF  Hid  THEN 
44R0  RFFP 

44R5  INPUT  E  . ' FR  X-I.AREI  * 

4490  RFFP 

4495  INPUT  “FNTFR  f  -I  APFI  “  .  vlab«»  !  S 

4500  pHO  IF 

4505  IP  Ono  2  THFN 

4510  PRINT  '9R  I , 2 i PU  PA  40  -14i* 

4515  PRINT  'IfliTiPR  -1.6.1  III  PR  1.7,0  PUsPP  .5,-4tl  BwniPR  .  S  ,  1  « * 

45^0  PRINI  "in-TirR  .5  -lilOdiFR 
4575  |F  l>m'0  THFN 
45 30  PPINT  -t.pl  IK)" 

4535  FI  -F 

4540  FP  TNT  "1  P  >  <  f  »" 

4545  FNO  ’F 

4550  FMP  IF 

4555  [F  Mpn  -=2  THEN 

45R0  IT  lun-0  THEN 

4^05  PRINT  "PP  I  571011  PA  40  »4||Pn  lU*ntf;R  I  I  .  5  j  PR  0.5  MlR7l5R  I.5,7,PR 

0.5  .  f  tl  PI' 

4570  FI  -F 

4r,7R  PRINT  "  5R  1.5.2IPU  PA  .14  .  - 1 4  tl  Ro  (Rln/hriPR  .5,.5ilfl.iPR  .5,-.5i" 
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4SP*  PRINT  "LBftiPR  .  5 , 1  i  SR  1  ,  I  .5tl.02 » SR  l.5,2iPR  .Sr-liLfl>i" 

4585  ENO  IF 

4590  END  IF 

4595  IF  Opo-0  THEN 

4fi00  IF  Iun-0  THEN 

4605  PRINT  “5R  1.5,2iPU  PA  -l2,40iOt  0,l*LBo  <W/«|PR  -I ,0.5iSR  I.I.StLBZtSff  1 
.  5 ,2 iPR  I  ,.5iLB >* 

4610  ELSE 

4615  PRINT  *SR  1 . 5 ,2 iPU  PA  -I2.32*DI  0,liL8q  (Btu/hriPR  - . S . . 5iL0 . iPR  .5. .Si' 
4620  PRINT  -LSftiSR  I.I.StPR  -1,.5iLB2iPR  1,.S»SR  T.S,2ilB>" 

4625  ENO  IF 
4630  END  IF 
4635  IF  Opo  0  THEM 
4640  IF  Iun-0  THEN 

4645  PRINT  "SR  1.5,2iPU  PA  -12 .38*01  0,IiL0h  <W/«iPR  -I,.SiSR  I,I.5iLB2iSR  I. 
5  , 2 i PR  .5 , .5i " 

4650  PRINT  "1.8.  i PR  .S.0iieK>" 

4655  ELSE 

4660  PRINT  "SR  1.5,2'PU  PA  -I2,20iOI  0,liL0h  (Btu/hriPR  - . 5 . . SiLB . iPR  .5. .Si" 

4665  PRINT  "LPftiPR  -I..S1SR  I.I.SiLBZiSR  l.5.2iPR  .S.-SiLB.iPR  .S(.SiLBF>* 

4670  ENO  IE 

4675  ENO  IF 

4680  IF  !  rll  *0  THEN 

4685  PRINT  "SR  1.5,2iPU  PA  50,-16  CP"  1  -I  EN(  F  lebe  If  )/2 1  "0  ilB"  1  Xlibe  1  ft " " 

4690  PRINT  'PA  -14.50  CP  0 . " i -1  EN( Ylabe I  * >/2«5 '6 1 *01  0  .  I iLB" 1 Vlabelti " " 

4695  PRINT  *CP  0.0  01* 

4700  ENO  IF 
4705  Ion-0 

4710  Reoeat :  • 

4715  411-1. E*G 

4720  Xul*-t.E+6 
4725  Icti'O 
4730  0EEP 

4735  INPUT  ’WANT  10  PLOT  OATA  FROM  A  FILE  < t-Y  ,0-N  >7"  ,04 
4740  IE  Olr-J  THEN 
4745  BEEP 

4750  INPUT  -ENTER  1  HE  NAME  OF  THE  DATA  FILE”  .0.1 1  left 
4755  ASSIGN  fFile  TO  O.filof 
4760  BEEP 
4765  8EFP 

4770  INPUT  -ENTER  THE  BEGINNING  RUN  NUMBER" ,Md 
4775  BEEP 

4700  INPUT  "ENTER  THE  NUMBER  OF  X-V  PAIRS  STORED" .Npairs 
4785  BEEP 

4790  INPUT  ‘CONNECT  DATA  WITH  LINE  »  t-V ,0*N  )7"  ,  Ic  I 

4795  BEEP 

4800  PRINTER  IS  I 

4805  PRINT  USING  ’4X, "-Select  *  symbol:'"" 

4010  PRINT  USING  *6Xf""l  St  nr  2  Plus  sign""" 

4815  PRINT  USING  -6¥r’"3  Circle  4  Sauare*"" 

4020  PRINT  USING  'GX.-’S  RqnMuS"*" 

4825  PRINT  USING  “6Xt""6  Right -aide-uo  triangle’"’ 

4030  PRINT  USING  "GX,-"7  Uo-s  1  de-down  irisnqle""" 

4835  INPUT  Syn 
4840  PRINTER  IS  705 
4845  PRINT  "PU  01" 

1850  IF  S*n-I  THEN  PRINT  "SM«* 

4055  IF  Svn-2  THIN  PRINT  "?«♦" 

4860  IF  Svn-3  THEN  PRINT  ’SMo’ 

4865  IF  Md  •  I  THEN 
4070  FOR  I’l  TO  f Md-I  ) 

4875  ENTER  0FileiYa,¥a 

4980  NEXT  I 

4P85  ENO  IF 

4890  FOR  1*1  TO  Noairs 

4895  ENTER  #Fl lei Ya .Xa 
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4900 
4  90S 
4910 
49 1 S 
4  920 
4925 
4930 
4935 
4940 
4945 
4950 
49SS 
4960 
4  965 
4970 
4975 
4  980 
4985 
4990 

4  995 
5000 
5005 
5010 
5015 
5020 
5025 
5020 
S075 
5040 
5045 
5050 
5055 
6060 
5065 
5070 1 
5075 
S060 
5085 
5090 
S095 
5100 
5105 
51 10 
51  15 
5120 
5125 
5130 
5125 
5M0 

5  US 
5150 
5155 
5160 
5165 
5170 
51 75 
5190 
5185 
5190 
5195 
5200 
5205 
5210 
5215 
5220 
5225 


IF  0po-»  THEN  Ya-Ya/Xa 

IF  0po-2  1HEN 

Q-Ya 

Ya-Ya/Xa 
Xa-fl 
END  IF 

IF  Xa<X  1 1  THEN  XU-Xa 
IF  Xa>Xul  THEN  Xul-Xa 
IF  Iun-1  THEN 
IF  0oo< 2  THEN  Xa-Xa*1.8 
IF  Ooo;0  THEN  Ya-Ya-.I76I 
IF  Opo-0  THEN  Ya-Ya« -2(7 
IF  0oo-2  THEN  Xa-Xa*.3l7 
END  IF 

X-L6T<  Xa/Xnin)«Sf.< 

Y-l6T(Ya/Yntn)»Sfv 
IF  3y*'3  THEN  PRINT  “SM" 

IF  Sy«'4  THEN  PRINT  ‘SR  1. 4.2.4' 

IF  lcl-0  THEN 
pniNT  “PA" ,X , Y , “ * 

ELSE 

PRINT  'PA' .X , i  , 'P0* 

END  IF 

IF  Sy*'3  THEN  PRINT  'SR  I .2,1.6" 

IF  Svn-4  THEN  PRINT  'UC2 .4  ,99  ,0 , -8 .-4 ,0 ,0 ,8 .4 ,0i '  I 

IF  Syn-S  THEN  FRINT  "UC3  .«  .99.-3  .-6  .-3,6  ,3  .6  ,3.-6i" 

IF  Syn-5  THEN  PRINT  'UC0  .5. 3  ,99 ,3 .-8 ,-6 .0 . 3 .0* ' 

IF  Svn-7  THEN  PRINT  MC0  . -5 . 3 ,99  .-3 .8  ,6  ,0  .-3  ,-8i ' 

NEXT  I 
PRINT  “PU* 

BEEP 

ASSIGN  OF.le  TO  • 

XII-V1J/I.2 
Xul-Xol • I .2 
6OT0  8040 
END  IF 

PRINT  -pu  St r 
BEEP 

INPUT  'WANT  TO  PLOT  A  POLYNOMIAL  M-Y  ,0-N  )»*  .Go.on 

IF  Po_on-1  THEN 

BEEP 

PRINTER  IS  I 

FRINT  USING  '4X. "Select  Ima  t*pe:'“' 

PRINT  USING  '6X  "0  5^ I  i d  line"' 

PRINT  USING  -RX,"I  OaeheiJ"  “ " 

PRINT  USING  “GX ,S  Longer  line  -  dash**’ 

INPUT  Ion 
PRINTER  IS  70S 
BEEP 

INPUT  'SELECT  <  0*L  IN  .  I -LOG<  X  ,  Y  ) )'  .  1 1  eg 
lornt- I 
CALL  Pclv 

FOR  r.-0  TO  C*  STEP  CsiCM 
Xa-X«m«  10  x* 

IF  xa  X|  J  OR  Xa  »Xul  THEN  5300 

frn- ten* 1 

Pu-0 

IE  Inn-1  THEN  tdf-I:n  MOO  2 
IT  Ion-2  THEN  Irtf-Icn  HOP  4 
IF  Ion-3  THEN  Mf-fcn  tinp  H 

IF  Ion-4  THEN  Mf-Icn  MUp  |6 

IF  fon*S  THEN  Mf-Icn  TWO  2P 

IF  Mf-l  IIIFN  Pm- I 

IF  Ooc-0  THEN  Ya-f  NPoIv<Xa  ) 

IF  Ooo-2  AMO  Ilog-0  THEN  > a-X.i /ENPn 1 y ( Xa  1 
IF  Opo-2  AND  I  log-1  IHFH  i  a-FNPo  1  v<  X«  I 
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5230  IF  Ooo-I  THEN  Ya-ENPoI.IXal-Xp 
S235  IF  Ya'Ymn  THEN  5x0a 
5240  IF  lun-  I  THEN 

5245  IF  Odo  2  THEN  Xa-xa-i.O 

5250  IF  Ooo>0  THEN  Ya*Ya«.l76l 

S2SS  IF  Opo-0  THEN  Ya-Ya-,317 

5260  IF  Opo-2  THEN  Xa->a«.3l7 

5265  END  IF 

5270  Y-LGTT  Ya/Ynin  >•$  f  v 

S275  X-L6TT  Xa.'Xni  n  )«STs 

5280  IF  Y:0  THEN  Y-0 

S28S  IF  YM00  THEN  GOTO  5300 

S290  IF  Pu-0  THEN  PRINT  "FA" .X , Y , "P0‘ 

5295  IF  Pu-I  THEN  PRINT  "PA"  Y  Y  -Pu‘ 

S300  NEXT  X., 

5305  PRINT  -PIT- 
S3 10  EHO  IF 
S3  IS  BEEP 

5320  INPUT  ‘ UAH I  TO  OUIT  <  | -Y ,0-N .  Iqt 

S32S  IF  rot  - 1  THEN  5335 

5330  GOTO  4715 

5335  PRINT  "PU  PA  0  0  SP0" 

5340  SUBEND 

5345  OEF  FNHstioot  hi  X  ,Boo  ,  1  sat  } 

5350  OITT  A(  51 .81 5  >  ,CT5>  ,DIS  ) 

5355  DATA  .20526, .25322. .319040, .55322, .79909  1.00256 
5360  OATA  . 74515, .72992 , .73199 , .71225 . .68472  .  .64197 
5365  OATA  41092 ,. 17726  .25142 , .54806 . .81916 , I .0845 
5370  OATA  71403, .72913. .72565, .696691 ,.665867, .61809 
5375  RE  AO  AT «  I  ,0(  •  I  ,C<  •  )  ,0T  •  ) 

5380  IF  Boo‘6  THEN  1-0, )0 
5395  IF  Boo-6  THEN  |>< 

5390  IF  Bop-10  THEM  1-5 
5335  IF  I  sat  -  I  THEN 
5400  Ho-EXP(A( I >»0< I  ).LOG< X  )l 
5405  ELSE 

5410  H»-EXP(C( I >*0< I >»LOG( X M 

S41S  END  IF 

5420  RETURN  Hs 

542S  FNENO 

5430  OEF  FNPoly(X) 

5435  CON  /Coly/  AT  I0.l0t.ci  10), 0(4  1  Nop , lornt  ,0po  Moo 
5440  Xl-X 
5445  Pol y -8T  0  > 

5450  TOR  I-l  TO  Noo 

5455  IF  I  log-1  THEN  XI-L06IX) 

5460  Poly«Poly«0T I )«X1‘I 
5465  NEXT  I 

5470  IP  Hog-1  THEN  Pol y-EXPTPoly T 
5475  RETURN  Poly 
5480  FMGMO 
SI8S  SUB  Poly 

5490  OIN  R(  10), ST  10  >,Sv(  12  1,  S- >  12  > ,x> ( l80>,Yv<  100  I 
S495  CON  'Coly/  AT  10, 10)  ,C(  101 .0(4  >  N.lprnt  ,0oo  .Hog 
5500  CON  /x-yy/  Xo( 5  I ,Yo<  5  > 

5505  FOR  i-o  TO  4 
5510  BID-0 
S5I5  NEXT  I 
5520  BEEP 

5525  INPUT  ’SELECT  (0-FILE  I • KEYBOARD , 2-PR06RAN  )•  .  In 
5530  In*  [pH  | 

5535  BtEP 

5r7<0  INPUT  -EM1W  NUMBER  OF  X-  1  PMFS"  No 
5515  IF  In- 1  THEN 
5550  BEEP 

5555  INPUT  “EM TER  DATA  FILE  NAME  *  0  filet 
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5560  BEEP 

5565  INPUT  ’LIKE  TO  EXCLUDE  OAJA  PAIRS  (I -Y  ,0-N >V , led 
5570  IF  ! ed» 1  THFN 
5S7S  BEEP 

55ee  INPUT  -ENTER  NUNGER  OF  PAIRS  to  BE  EXCLUOEO' .Ipex 
5585  END  IF 

5590  ASSI6N  *Ftle  TO  0_ftle» 

5595  ELSE 
5S00  BEEP 

5605  INPUT  *UANT  TO  CREATE  A  DATA  FILE  < 1 -V r0-N >?-  .Yes 
5610  IF  Yes-1  THEN 
5615  BEEP 

56:0  INPUT  ‘GIVE  A  NAME  FOR  DATA  FlLE*,0_fitef 
5675  CREATE  BOAT  O.file*  5 
5630  ASSIGN  0fiie  TO  OJlle* 

5635  END  IF 
"640  END  IF 
545  BEEP 

5650  INPUT  *ENTFR  THE  ORDER  OF  POL  LNONIAL  * ,M 

5655  FOR  1-0  TO  N*2 

5660  Sy(I>-0 

5665  Sa(  I  >-0 

5570  NEXT  I 

5675  IF  led- I  AND  I«-1  THEN 

5630  FOR  I- I  TO  toex 

56F5  ENTER  0File»X,Y 

5690  NEXT  I 

5695  ENP  IF 

5700  FOR  I- I  TO  Mo 

5705  IF  Tn-f  IHEN 

5710  IF  Ooo-Z  THEM  ENTER  CFiie«x,V 

5715  IF  Ooo'Z  THEN  ENTER  0Fiie<Y,X 

5730  IF  I  log- 1  THEN 

5735  xt-x/Y 

5730  X-LOG(X) 

5735  Y-LOGI X I ) 

5710  END  IF 
5745  END  IF 
5750  IF  I*-2  THEN 
5755  BEEP 

5760  INPUT  -ENTER  NEXT  X  ■>’  PAIR*,X.Y 

5765  IF  Yes- 1  I  HEM  OUTPUT  0File«X,Y 

5770  6 NO  IF 

5775  IF  I*  3  THEN 

5780  X..(tJ-X 

s?es  Yv<I)-Y 

5790  ELSE 

5795  X-Xod-M 

5300  Y-YoII-n 

5BP5  END  IF 

5910  R(0)-Y 

SB  1 5  Sy<0>-Sv<0)*Y 

5830  S< 1  >-X 

5935  S.*T  I  )-ShI  I  )♦« 

5830  FOR  I- I  TO  N 
5835  R<J  HdlJ-HM 
5910  9»(J>-Sy(J)*R(J) 

5045  ME X  T  J 
5850  FOR  J-Z  10  N*2 
5555  Si J  >*S<  J- I  i *X 
5060  5  .  <  j  w s ^  <  J  >»«■><.!> 

5065  ME xj  J 
5070  NEXT  1 

SB75  IF  Yes-1  AND  I«-2  THEN 
5080  BEEP 

5995  PRINT  USING  "I2X.D0,""  x-Y  oairs  were  stored  in  file  ,  I0A* iNo  0_f 1 le* 
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5890  END  IF 
5895  S.-.<0>-No 

5900  FOR  1-0  TO  N 
5905  C< I >-Sy i I  ) 

5910  FOR  J-0  TO  N 
59  IS 

5920  NEXT  J 
5925  NEXT  1 
5930  FOR  1-0  TO  N-l 
5935  CALL  Divide* l ) 

5940  CALL  Subtract! !♦!  ) 

5345  NEXT  I 

5950  B<M>-C<NKA<N,N> 

5955  FOR  1*0  TO  N-l 
S960  8<  N- 1 -I )»Cf  N-l-I > 

5955  FOR  J-0  TO  1 

S970  0< N- I -I  )-Bf M- I -I >-A<N-1 -I ,N  J  )»B( N- J  J 

S975  NEXT  J 

5300  B<N-I-1  >*B<N-1-I  >/A<N*  1-!  M-1-l  ) 

S99S  NEXT  I 
5990  'PRINTER  IS  701 
5995  *  PRINT  8< • ) 

6000  'PRINTER  15  70S 
5005  IF  Iornt-0  THEN 

6010  PRINT  USING  *  J  2X  ,  * "EXPONENT  COFFF ICIENT * " ’ 

6015  FOR  1*0  TO  N 

6020  FRrNT  USING  H I SX  ,00  .5X  ,flD .  70E  “  i  I  ,B(  I  ) 

6025  NEXT  I 
6030  PRINT  ‘  - 

6335  PRINT  USING  •  I2X  .‘‘"DATA  POINT  X  Y  V(  CAL  Cl  H.  ATEO  i  OISCR 

EF  ANCY" 

6040  FOR  I- I  TO  No 
6045  Yc-8<  0  > 

6050  FOR  J-l  TO  N 
6055  Yc- Yc  »B<  J  > • X* ( I  )  '  J 
6060  NEXT  J 

6065  D-Xyll)  Ye 

6070  PRINT  USING  *  ISX  ,30  ,4X  .  4( MD50E  .  I X  >"  1 1  ,X.\<  I  > ,  Y>(  I  > .  Yc  .0 

6075  ME x f  I 

6090  END  IF 

6005  ASSIGN  tfFile  TO  • 

6090  SUBFNO 
6095  SUB  Dl  vide! M  > 

6100  CON  /Coly/  A(  10, 10  )  ,C<  10)  ,0<4  )  ,N,Iprnt  ,0do  .IIoo 
6105  FOR  I-N  TO  N 
6110  Ao-A< 1 ,M> 

61  IS  FOR  J-M  TO  N 

6120  ATI  ,  J  >'A<  I  J  >  'Ao 

6125  NEXT  J 

6130  C<  I  i»C<  I  >/Ao 

6135  NEXT  I 

6U0  SUBEND 

6145  SUB  Subtree t <  F  ) 

CI50  CON  /Coly/  AM0, 10  )  .C(  10  >  Bl  4  )  ,N  .  Iprnt  ,0oo  ,  1  leg 

615S  FOR  1 ”K  TO  M 

6160  FOR  J-K-l  TO  N 

6165  A«  I  .  J  )»A< K *  I  J  >-A<  I  ,  J  ) 

6170  NEXT  J 

6175  C<  I  )-C(K  I  )  C(  t  ) 

6190  NEXT  I 
6185  SUBFNO 
6  I  9*3  SUB  PI  in 

6195  CON  /Colv/  A(  !0,I0>,CM0)  B(4  >  N.Iornt  .Ooo.I  log 
6200  COM  /X,vy/  X.<S».Y,(5) 

6206  PRINTER  IS  70S 
6210  BEEP 
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6215  INPUT  ‘SELECT  <0-h/h0*  aawe  lube  .  I  "ht  HF  >/h(  an  1‘ ,  Irt 
6220  BEEP 

6225  INPUT  ‘WHICH  Teat  ( I -6.7  ,0--2 . 2 >" . leal 

6220  Xntn“0 

6235  Xme>.-I0 

6240  Xatep-2 

G24S  IF  Irt-0  THEN 

6250  Vmn-0 

6255  Yima-1.4 

6260  Yaleo-.2 

6265  ELSE 

6270  Ymn-0 

6275  Y«a.-.-IS 

6290  Xateo*5 

6295  END  IF 

6290  BEEP 

629S  PRINT  ‘INiSPliIP  2300 ,2200 .8300 ,68001 ‘ 

6300  PRINT  ‘SC  0 , 100 ,0 , 100t  TL  2  Oi* 

6305  Xnaa-Xnm  > 

6310  Sfy» 100/1 Yna--Yntn ) 

6315  PRINT  ‘ PU  0.0  PO- 

6320  FOR  Xa*Xnin  TO  Xna.<  ST EP  Xateo 

632S  X»(Xa-Xntn)>$f« 

6320  PRINI  ‘PA‘ i X , ‘ .0i  XT i ‘ 

6335  NEXT  Xa 

6340  PRINT  ‘PA  100  ,0iPUi ‘ 

6345  PRINT  PU  Pfl  0 .0  P0‘ 

6353  FOR  Xa*Ynln  TO  Xna*  STEP  Yaleo 
6355  Y-<Ye-Ymni»Sfy 
6360  PRINT  ‘PA  0  ‘ |Y ,‘YT‘ 

6365  NEXT  Yj 

6370  PRINT  ‘PA  0,100  fL  0  2‘ 

6375  FOR  Xa-Xmn  TO  Xna.-,  STEP  Xateo 
6380  X><  xa-Xmn  >*ST  < 

6395  PRINT  ‘PA‘ iX  ,  1001  XT* 

6330  NEXT  Xa 

6395  PRINT  ‘PA  100,100  PU  PA  100  0  PO" 

6400  FOR  la'Vxin  TO  Viaa-  STEP  Yateo 

6405  Y-( xa-Ywib  i-Sfy 

6410  PRINT  ‘PO  PA  100.- ,Y,‘YT‘ 

6416  NEXT  xa 

6420  PRINT  ‘PA  100,100  FU‘ 

6425  PRINT  ‘PA  0,-2  SR  1 .5,2‘ 

6430  FOR  Xa-Xmn  TO  Xna*  STEP  Xateo 
6435  X-H  xa-Xnin  >»Sf  > 

6440  PRINT  ‘PA*lX  ,‘,0i‘ 

6445  PRINT  ‘CP  -2  ,-l iU8‘ iXai  “ 

64S0  NEXT  xa 

6455  PRINT  "PU  PA  0,0‘ 

6460  FOR  Xa-XAln  10  Yno*  STEP  Xateo 
6465  IF  A0S(  Ta)  I.E-S  THEN  x.a»0 
6470  Y*1 Ya-Xnln )*Sfy 

6475  PRINT  -PA  0,‘iX,‘* 

6480  PRINT  ‘CP  -4  ,-.2SiLB‘ iXai "• 

6465  NEXT  Y.J 

6490  X|abel»*‘0tl  Percent" 

6495  IF  lrt*0  THEN 
6500  Xlabel*-  A  b0V 
6505  ELSE 

6510  xlabei«*‘NThaetoolh‘ 

6515  END  IF 

6520  PRINT  "SR  1.5.2  (PU  PA  50  -10  CP"  t  EH*  X  label  *  1/2i'0iI8’iX  labe  1  •  i 

6525  PRINT  ‘PA  -11,50  CP  0  ‘  i -l  EII(  v  labe  I*  >' 2 -S-Si  "01  0 , 1  il  B‘  i  x  label «;  “ 

6530  PRINT  'CP  0  0‘ 

653S  lon-0 
6540  8FEP 
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6545  INPUT  'WANT  TO  PLOT  DATA  FROM  A  FILE  <  I -V  0-N  >?*  .01- o 

6550  Icn-0 

6555  IF  mg-)  THEN 

65G0  BEEP 

G565  INPUT  ‘ENTER  THE  NAhE  OF  THE  DATA  FILE* .0  filet 
6570  SEEP 

6575  INPUT  ‘SELECT  (0-L  INEAR.  1  -LOG<  X  ,Y  >"  .  Hog 
6580  ASSIGN  tFile  TO  O.flle* 

65B5  BEEP 

6590  INPUT  ‘ENTER  THE  BEGINNING  RUN  NUMBER* ,Md 
6595  BEEP 

6600  INPUT  'ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED*  .NpalPs 
6605  BEEP 

6610  INPUT  ’ENTER  DESIRED  HEAT  FLUX*  ,0 

6SIS  BEEP 

6620  PRINTER  IS  I 

6625  PRINT  USIN6  ‘4X, '"Select  a  synbol ! "" " 

6630  PRINT  USING  ‘4X,""|  Star  2  Pius  siQn‘"" 

6635  PRINT  USING  *4X,""3  Circle  4  Square*** 

6640  PRINT  USING  ‘4X,**5  Ro«bus**“ 

6645  PRINT  USING  *4X.‘*6  Right -side-ug  triangle""" 

6650  PRINT  USING  *4X,*"7  Up-slde-down  triangle*** 

665S  INPUT  Svii 
6660  PRINTER  IS  705 
6665  PRINT  *PU  DI" 

5670  IF  Syn-1  THEN  PRINT  ’SM*“ 

6675  IF  Sv«-2  THEN  PRINT  "SM»" 

6GeO  IF  Sy«-3  THEN  PRINT  ‘SMo* 

8665  Nn»4 

6690  IF  I  log-1  THEN  Nn-I 
6695  IF  Md.'  1  THEN 
6700  FOR  I-l  TO  (Md-t ) 

6705  ENTER  tFiletXa.Ya 
6710  NEXT  I 
6715  END  IF 
6720  0 1 *0 

6725  IF  1  log-1  THEN  0-1.06(0  > 

6710  FOR  I-I  TO  Noairs 
6735  ENTER  SFileiXa  ,0< •  ) 

6740  Ya«fl<  0 ) 

6745  FOR  K » I  TO  tin 
67S0  Ya-Ya+0(k  )«0‘K 
6755  NEX  T  K 

67G0  IF  I  log- 1  THEM  Ya-EXP(Va) 

6765  IF  I  log-0  THEN  Ya-QI/Ya 

5770  IF  Irt-0  I  HEN 

6775  IF  Xa-0  I HEN 

67B0  Yo-Ya 

6795  Ya-I 

6790  ELSE 

6795  Ya- Ya/Yo 

6600  ENO  IF 

6805  ELSE 

6810  Hsn-FNHsnoot h(  0 fXa , Isat ) 

6015  Ya-Ya/Hsn 
6020  ENO  IF 
6025  in  I-l  )-Xa 
6030  Yyl I-l  )-Ya 
6835  X*( Xa-Xnin  >«Sf « 

G040  Y  »<  Ya-Ym  n  )  ■  S  f  y 

6045  IF  Svn  •  3  THEN  PRINT  *5M‘ 

6050  IF  5v«  4  TIICN  FPINT  “SR  1.4, 2. 4* 

6055  PRINT  *PA ‘  .X  .Y  .  *  * 

6060  IF  Syn-3  THEN  PRINT  *  SR  1. 2.1.6" 

6665  IF  Syn-4  THEN  PRINT  "UC2.4.99  0  . -0  .  -4  .0  ,0 , 9  4  . 0  .  i ‘ 

6870  IF  5  vci-5  MIEN  PRINT  "Uf‘3.0  99  .  -  3  . -6  .  -  3  6  . 3  .6  . 3  . -fi  i  ‘ 
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S87S  IF  Sy«-6  THEN  PRINT  *UC0 ,5. 3 , 99 ,3 . -0 , -6 . 0 , 3 , 8 r * 

6880  IF  Syn-7  THEN  PRINT  *UC0  ,-5. 3, 99, -3,8. B, 0,-3, -8 1" 

6085  NEXT  I 
6090  BEEP 

6095  ASSIGN  8Ftle  TO  • 

6900  END  IF 
6905  PRINT  *PU  $N’ 

6910  BEEP 

691S  INPUT  ‘WANT  TO  PLOT  A  POLYNOMIAL  <1 -Y , 0-N >7 ‘ ,0V p 
6920  IF  Oko-t  THEN 
6925  BEEP 

6930  INPUT  ‘SELECT  (  0-LINEAR , l-LOGl X  ,Y  )  >“  .  I  log 
693S  lornt-l 
6940  CALL  Poly 

6945  FOR  Xa-Xnin  TO  Xnax  STEP  Xateo/25 
6950  Icn*Icn*l 
6955  Ya-FNPoiy{ Xa ) 

6960  Y-<  Ya-Ymn  ) •  S  f y 

6965  X-<Xa-Xmn>«Sfx 
6970  IF  v.,0  THEM  Y-0 

5975  ir  Y  100  THEN  GOTO  7025  # 

6980  Pu-0 

6985  IF  Ion-1  THEN  Idf-Icn  MOO  2 

6990  IF  Ion-2  THEN  Idf-Icn  MOD  4  / 

6995  IF  Ion-3  THEN  Idf-Icn  MOO  8 

7000  IF  Ion-4  THEN  Idf-Icn  MOO  16 

7005  IF  Ion-5  THEN  Idf-Icn  MOD  32 

7010  IF  Idf-l  THEN  Pu-f 

7015  IF  Pu-0  THEN  PRINT  *PA"  ,X  f,*PO’ 

7020  IF  Pu-f  THEN  PRINT  ‘PA "  .  X  ,  Y  .  “PU‘ 

7025  NEXT  Xa 
7030  PRINT  *PU‘ 

7035  Ipn-Ipn+I 
7040  GOTO  6540 
7045  ENO  IF 
7050  BEEP 

7055  INPUT  ‘WANT  TO  OUIT  < I -Y  ,0-N >7*  ,  Iqui l 
7060  IF  Iquit-l  THEN  7070 
7065  GOTO  6540 
7070  PRINT  ‘PU  SP0‘ 

7075  SUCEND 
7080  SUB  Stats 
7065  PRINTER  IS  701 
7090  J-0 

7095  K-0 

7100  BEEP 

7105  INPUT  ‘PLOT  FILE  TO  ANALYZE’’  ,F l let 
71 10  ASSIGN  OF i 1 e  TO  File* 

71)5  BEEP 

7120  INPUT  ‘LAST  RUN  No ’( 0-9U1T  >*  (Nn 

7125  IF  Nn-0  THEN  7305 

7130  Nn-fln-J 

7135  S*-0 

7140  Sy-0 

7145  S:-0 

7150  5*9-0 

7155  Svs-0 

7160  Sc  s-0 

7165  FOR  I-l  TO  Nn 

7170  J-JM 

7175  ENTER  OF : l e : 0 , T 

7180  H-O/r 

7185  S*«S.*»0 

7190  S • s-S ■ s  *0  '2 

7195  Sy-Sy*T 

7200  Sy5*Svs»T-2 


.  ->  •  .  ■ 
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Sz-SztH 

Si»-Si»tH*2 

NEXT  I 

Oava-Sx/Nn 

Tava-Sy/Nn 

Hava-Si/Nn 

Sdavq-S0RIA8S<  INn'5-J-S.  2  >/(Nn*(Nn-l  ) ) )) 

Sdavt-S0R(A8S<  <Nn*Sy»-Sy  2  )/INii*(Nn-l  ) t ) * 

Sdevh-SOflt ABSI (Nn*5zs-Sz  2  >/<Nn*TNn-t  >)>) 

Sh-l00*S da  vh. -’Ha  ye 

Sa-l00*Sde-'q/Qave 

St-i00*SdevL/Tav« 

IF  K-l  THFN  7295 
PRINT 

PRINT  USING  " 1 1 X . ‘ “DATA  FILE: "* . I4A" iFi |qt 
PRINT 

PRINT  USING  "I  IX," “RUN  Htube  SdevH  0 


SdevQ  Thatab  SdevT" 


7290  K-l 

7295  PRINT  USING  *  I  I X  ,DO  ,2< 2X ,0 . 30E  ,  I X .30. 20 > ,2X  ,00 . 30 , 1  X  ,3D . 2D" i J ,Hav« ,Sh .Gave 
,Sq ,Tave ,St 
7900  GOTO  71  IS 
7305  ASSIGN  OFilel  TO  • 

7310  PRINTER  IS  I 
7315  SUBEriO 
7320  SUB  Coef 

7325  COM  JCpiyS  AT  10, 10  )  ,CI  I0>,B(4  >,N,Iprnt  ,0oo, U09 
7330  PCEP 

7335  INPUT  "GIVE  A  NAME  FOR  CROSS-PLOT  FILf.CpfS 
7340  CREATE  BOAT  CoT*,2 
7345  ASSIGN  «Flla  TO  Cqf» 

7350  BEEP 

7355  INPUT  "SELECT  I  0-1.  INEAR ,  I  -I  OGI  X  ,X  I  )*  ,  1  log 
73G0  BEEP 

73S5  INPUT  ‘ENIER  OIL  PERCENT  ( -  I -STOP ) - ,Boo 

7370  IF  Bop  0  THEN  7390 

7375  CALL  Poly 

73GO  OUTPUT  0Ftl«t8op ,81 • > 

738S  GOTO  7360 

7390  ASSIGN  0Eile  TO  • 

7395  SUBENO 
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CO HPOTEB— CONTROLLED  VALVE  PROGRAM 

The  computer-controlled  valve  was  controlled  by  an 
Octagon  Systems  STS-2A  microcontroller.  The  SYS-2A  is  a 
complete  computer  system  on  a  single  card  and  requires  only 
a  5  volt  supply  for  operation.  The  SYS-2A  has  provisions  for 
4,  2.5  volt  A/D  data  inputs  and  8  high  current  digital 
outputs.  The  card  was  configured  as  shwon  in  Table  4  below. 


TABLE  4 

Configuration  of  SYS- 21  I/O  Channels 


Data  Input  Channel 
3#A10  *  8 
a#AlO  =  9 


Purpose 

Yalve  position 

Liquid  R-114  Temperature 


Data  Output 
a#A00 
2#aoi 


Channel 


Purpose 

Shut  Valve  Relay  (0=stop,  1=shut) 
Open  Valve  Relay  (0=stop,  1=open) 


The  proportional-integral-derivative  control  program 
below  was  written  in  NSC  "tiny  BASIC."  The  SYS-2A  microcont¬ 
roller  operator's  manual  includes  an  appendix  of  tiny  EASIC 
commands.  Remarks  are  provided  below  that  were  not  in  the 
actual  program  (to  minimize  execution  time)  for  clarity.  The 
program  was  stored  in  EEPROH  added  to  the  card  after  program 
debugging  was  complete. 

10  Let  A=15;Let  B=5: Let  C=10:Let  F=0:Let  1=0  Initialize 
15  Let  3tA00=0:Let  3#A01=0  Stop  Valve 

20  Print"Set  Constants (1=Y, 0=N) ?":Input  Q 
30  If  (Q=0)  then  GOTO  60 
40  Prim t"A#B,C,G?": Input  H,J,K,L 
50  Let  A=H:Let  B=J:Let  C=K: Let  G=L 
60  Print"Input  Desired  Temp (-1=End) ":Input  R 


Enter  Temp 


70  Let  3#A10=9 : Delay  2:Let  M=16*3#A1 1+9*412/16  Head  Channel  9 
80  Let  a#Al0=8:Delay  2:Let  P=16*3#A1 1+3#A12/16  Bead  Channel  8 
90  Print"Temp=,,#S,,,Pos=,,,P 

100  If(B<0)then  GOTO  800  Exit  Program 

200  Let  3#A10=8:Delay  2:Let  P=16*a#A11+a#Al2/16  Read  Channel  8 
210  Let  S=0  Are  10  Readings 

220  For  N=1  to  10  Step  1 

230  Let  3#AlO=9:Delay  2:Let  M=16*3#A11+a#Al2/1 6  Read  Temp 
240  Let  E=M-R: let  S=S+E  Error ; average 

250  Next  N 


260  let  E=S:Let  D=S-F:let  F=E 

262  If  (E<0)  then  1=1-1 

264  If  (E>0)  then  1=1+1 

266  If  (I <-31 000)  then  I=-31000 

268  If (I>31000)  then  1=31000 

270  Print"N=",H#"E="#I,"D=",D,  "I=",I 

300  If  (E<- 1000)  then  GOTO  600 

310  If  (E> 2000)  then  GOTO  700 

400  let  7=  (E/A)  ♦  (B*D)  ♦  (I/C) 

410  Print"V=",V,"ET=",E/A,"DT=",B*D,"IT=",I/C 
420  If  (V>0)  then  GOTO  500 
430  Let  V=V*(-1)*G 


Error; Derivative 
Integr al 
Integral 
Binary  Limit 
Binary  Limit 


Control  Band 
Control  Band 
PID  Value 
Debug  Tool 
Above  Desired 
Shut  Faster 
Chech 

Delay  Function  limited  1-1040  entry 
Slow  Valve  near  end  of  travel 


440  If  (V<1)  then  V=1 
450  If  (V>1 040)  then  V=0 
460  If  (P<150)then  V=15 
465  If  (P<50)  then  V=10 

470  If {P<20)then  GOTO  200  End  of  Travel  Limit 

480  Let  3#A00=1: Delay  (V)  :Let  a#AO0=O  Shut  Variable  Time 

490  GOTO  200 

500  If  (P<2500)  then  GOTO  200  End  of  Travel  Limit 
510  If  (V<1)  then  V=1  Check 

520  If  (V>1040)  then  V=0  Delay  Function  limited  1-1040  entry 
530  Let  3#A01  =  1 : Delay  (V) : Let  a#A01=0  Open  Variable  Time 

540  GOTO  200 

600  If  (P<20)then  GOTO  200  End  of  Travel  Limit 
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Oat  of  Bound  Warning 
Bax  valve  shutting  speed 
Slow  near  end  of  travel  Unit 


610  Print"600  loop" 

620  Let  7=0 

630  If  (P<250)then  7=50 
640  If  (P<150)  then  7=15 
645  If  (P<50) then  7=10 
650  let  3#A00*  Is  Delay  (7)  :  let  3#A00=Q  Shut  7alve  Fully 

660  let  a#A10=8 : Delay  2: Let  P=16*a#A1U3*Al2/16  Bead  Channel  8 
670  GOTO  600 

700  If (P>2500)  then  GOTO  200  End  of  Travel  limit 

710  Print"700  loop"  Out  of  Bound  Warning 

720  let  a#A01= 1 : Delay  0:let  a#A01=1  Open  Valve  Fully 

730  let  3#A10*8:Delay  2: let  P= 16»a#A11+a#A12/16  Bead  Channel  8 

740  GOTO *700 

800  END 


1PP EH Dll  D 

EXMBLES  OF  BEPBESEET1TIYE  D1TX  HOMS 

The  tvo  printouts  belov  are  sanples  of  the  output  of 
DBP2.  The  first  printout  (WH79)  is  for  the  saooth  tube  at 
high  heat  flux  with  0  percent  oil.  The  second  printout 
(HF101)  is  for  the  High  Flux  tube  under  siniliar  conditions. 


Month,  data  and  tine  :03:2S: 10:02: 18 

NOTE:  Progran  nane  :  QRP2 
Olsk  nuntaer  -  02 
Old  file  nane:  WH79 

This  data  set  taken  on  :  02:10:12:50:15 
Tube  Type:  4 

Data  Sat  Nunber  *  I  Sulk  Oil  I  -  0.0  02:10:13:16:10 

TC  No:  I  2  3  4  5  6  7  8 

Teno  :  16.42  14.69  14.61  15.00  14.67  14.91  16.79  14.74 

Tua  Tltod  Tliad2  Tvaor  Psat  Tsuno 

IS. 23  -2.26  -2.20  -.63  -6.31  -17.7 

Thetab  Htube  Odo 

I6.93S  S.704E«03  9.6606*04 

Oata  Set  Nunber  •  2  Bulk  OH  X  •  0.0  02:10:13:19:35 

TC  No:  I  2  3  *  5  6  7  0 

Teno  :  16.40  14.60  14.54  14.93  14.52  14.74  16.69  14.67 

Tua  Tliad  Tliad2  Tvaor  Psat  Tsuno 

IS. 15  -2.27  -2.10  -1.44  -6.32  -17.1 

Thetab  Htube  Odo 

16.866  5. 745E+03  9.6896*04 

Oata  Sat  Nunber  -  3  eulk  Oil  X  •  0.0  02:10:13:19:56 

TC  No:  I  2  3  4  5  6  7  S 

Teno  :  16.48  14.56  14.52  14. 9S  14.50  14.67  16.82  14.67 

Two  Tliqd  Tliod2  Tvaor  Psat  Tsuno 

IS. 12  -2.27  -2.13  -I .SI  -6.32  -17.0 

Thetab  Htube  Odo 

16.840  S.760E+03  9.7016*04 

NOTE:  03  X-Y  oalrs  were  stored  in  plot  data  file  PUH79 

DATA  F ILE :PUH79 

RUN  Htube  SdevH  Odo  SdevO  Thetab  SdevT 

3  5.736E+03  .51  9.683E*04  .22  16.880  .29 


Month,  date  and  tine  :03:2S: 10:03:26 

NOTE:  Program  name  :  0RP2 
Oisk  number  *  02 
Old  file  name:  HFI01 

This  data  set  taken  on  :  02: tS: 19:06:52 
Tube  Type:  5 

Data  Set  Number  -  I  Bulk  011  X  -  0.0  02:16:15:31:02 

TC  No:  I  2  3  4  5  6  7 

Temp  :  5.S2  2.63  3.76  4.57  3.93  4. 52  7.31  4 

Twa  Tllqd  Tllqd2  Tvapr  Psat  Tsunp 

4.55  -2.14  -1.87  .37  -6.20  -15.7 

Thetab  Htube  Qdo 

3.370  2.80IE+04  9.440E+04 

Oata  Set  Number  -  2  Bulk  Oil  X  -  0.0  02:16:19:31:16 

TC  No:  I  2  3  4.5  6  7 

Temp  :  5.54  2.65  3.76  4.57  3.93  4.53  7.31  4 

Twa  Tllqd  Tliqd2  Tvapr  Psat  Tsunp 

4.56  -2.14  -1.90  .32  -6.20  -15.6 

Thetab  Htube  Odp 

3.379  2.792E+04  9.433E+04 


Oata  Sat 

Nunber 

•  3 

Bulk  011  X  - 

0.0  02:16:19:31:27 

TC  No: 

1 

2 

3 

4 

5  6  7 

Teno  : 

s.ss 

2.63 

3.76 

4.57 

3.93  4.53  7.31  4 

Twa 

Tllqd 

T 1 lad2 

Tvapr 

Psat 

Tsuno 

4.56 

-2.14 

-1 .91 

.26 

-6.20 

-15.5 

Thetab  Htube  Qdo 

3.369  2. 0076+04  9.4566+04 

NOTE:  03  X-Y  pairs  were  stored  in  plot  data  file  PHFI01 
OATA  FILE:PHFI01 

RUN  Htube  SdevH  Qdo  SdevQ  Thetab  SdevT 

3  2 . 900E+04  .27  9 . 443E+04  .12  3.373  .16 


APPBHDI1  E 
0BCEBT1INTI  A Nil ISIS 


The  uncertainty  of  the  heat-transfer  coefficient  at  37 
kw/a2  and  5  ki/a2  of  runs  WH79  and  HP1Q1  are  analyzed  below. 
The  analysis  is  based  on  the  Kline-HcClintock  [Bef.  24] 
■ethod  of  uncertainty  analysis. 

The  heat-transfer  coefficient  is: 

(E.1) 

/ 

'  T  -  T  *  Twt  -  Qc  _  x  (E.2) 

wo  sat  2  ir  k  L  sat 

where 

h  *  heat-transfer  coefficient 
gc  *  heat  flux  corrected  for  end  losses 
^wo  3  average  outer  wall  teaperature 
Tsat  3  saturation  teaperature 

=  average  inner  wall  teaperature 
Qc  =  heat  input  corrected  for  end  losses 
D2  =  outer  wall  diaaeter  of  tube 
=  inner  wall  diaaeter  of  tube 
k  -  theraal  conductivity  of  wall 
1  *  length  of  heated  surface 

let 

_  Qc  In  (D2/D1) 

'i.H -  !E-3> 
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According  to  Kline  and  HcClintock,  the  uncertainty  in  the 
heat-transfer  coefficient  is: 


by  neglecting  the  ericr  from  the  logarithmic  tern,  because 
it  is  small  compared  to  the  other  terms,  the  uncertainty  of 
the  Fourier  term  (F)  can  be  estimated  as: 


and 

“  qc  »  d2  l  (E.6) 

which  has  an  uncertainty  of 


(E.7) 


Table  5  lists  the  various  terms  of  eguations  (E.1)  through 
(E. 7)  assuming  all  8  wall  thermocouples  are  used  to  get  the 
average  wall  temperature  ?  ^ .  The  data-reduct ion  pregram 
(DRP2)  used  this  method  to  calculate  the  heat-transfer 


1.T 


► 


TABLE  5 

Uncertainty  Analysis  Terns  Using  8  Thermocouples 

1 


„* 

File 

Heat  Flux 

WH79 

37  kW/a2 

WH79 

5  kfl/m2 

HF101 

37  kfl/m2 

HF101 

5  kH/n  2 

i 

6QC 

0.007 

0-  009 

0.006 

0.  007. 

6k 

"k 

0.15 

0.15 

0.33 

0.33 

l r» 

M 

5L 

T 

0.0005 

0.0005 

0.0005 

0.0005 

F  term  (°C) 

0.201 

0.0  24 

1.  316 

0.  170 

i 

S?  (OC) 

0.030 

0.004 

0.434 

0.056 

fwi  <°C> 

11.46 

5.63 

0.79 

-1.43 

a 

Tsat  <°C) 

-2.21 

-2.23 

-2.21 

-2.  18 

«Twl 

"^wo  “  "^sat 

0.032 

0.073 

0.290 

0.  126 

i 

6F 

^wo  “  ^sat 

0.002 

0.0005 

0.259 

0.097 

^aat 

^wo  ”  Tsat 

0.003 

0.008 

0.015 

0.  019 

i 

Sqc 

qc 

0.003 

0.007 

0.008 

0.003 

6h 

~h 

0.032 

0.074 

0.389 

0.  160 

i 

h  (W/a2  K) 

2660 

540 

22300 

8400 

» 


t 
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coefficient.  Table  6  lists  the  same  terms  assuming  only  the 
center  4  wall  thermocouples  (2,  3,  5,  and  6)  are  used  to 
calculate  the  heat -transfer  coefficient.  Comparing  Tables  S 
and  6  shoes  the  effect  of  the  axial  and  circumferential  wall 
temperature  distributions  on  the  uncertainty  of  the 
heat-transfer  coefficient. 

The  constraining  error  of  the  smooth  tube  is  the  uncer¬ 
tainty  in  the  wall  temperature.  Removing  the  effect  of  the 
axial  wall  temperatue  distribution  reduces  the  uncertainty 
by  1  percent.  The  large  wall  superheat  of  the  smooth  tube 
contributes  to  the  small  magnitude  of  the  uncertainty  terms. 
The  constraining  error  of  the  High  Flux  tube  is  also  the 
uncertainty  in  the  wall  temperature,  but  the  uncertainty  of 
the  thermal  conductivity  (part  of  the  F  term)  is  about  the 
same  magnitude  an  results  in  the  larger  overall  uncertainty 
of  the  High  Flux  data.  The  small  wall  superheats  of  the 
High  Flux  tube  also  amplify  the  magnitudes  of  uncertainty 
terms.  Removing  the  effect  of  the  axial  wall  temperature 
distribution  makes  the  uncertainty  in  the  thermal  conduc¬ 
tivity  of  the  copper- nickel  High  Flux  tube  the  constraining 
uncertainty.  The  axial  temperature  distribution  is  respon¬ 
sible  for  about  25  percent  of  the  wall  uncertainty  term,  but 
again  the  combined  effect  of  the  uncertainty  of  the  thermal 
conductivity,  and  low  wall  superheats,  does  not  make  it 
wholely  responsible  for  the  large  uncertainty  of  the  High 
Flux  data,  (lore  accurate  data  could  be  obtained  on  the  High 
Flux  surface  by  using  a  solid  copper  tube  without  a  large 
uncertainty  in  the  wall  resistance  or  an  axial  wall 
temperature  variation. 


TABLE  6 

Oncer tainty  Analysis  Terns  Osing  Center  4  Thermocouples 


File 

OH79 

SH79 

HF101 

HF101 

Heat  Flux 

37  kH/m2 

5  kW/m2 

37  kW/m2 

5  kW/m2 

6QC 

"q7 

0.007 

0.009 

0.006 

0.007 

6k 

T 

0.15 

0.15 

0.33 

0.33 

6L 

0.0005 

0.0005 

0.0005 

0.0005 

F  tern  (°C) 

0.201 

0.024 

1.316 

0.  170 

SF  (°C) 

0.030 

0.004 

0.434 

0.056 

f  .  (OC) 

11.34 

5.80 

0.44 

-1.49 

W1 

T  (°C) 

-2.21 

-2.23 

-2.21 

-2.  18 

sat 

«Twl 

0.015 

0.063 

0.213 

0.097 

^uo  "  ^sat 

6F 

0.002 

0.0004 

0.325 

0.108 

Two  ~  Tgaf 

®Tsat 

0.003 

0.008 

0.019 

0.021 

Two  “  Tgat 

6qc 

0.003 

0.007 

0.008 

0.003 

4c 

6h 

~ 

0.016 

0.064 

0.389 

0.  146 

h  (»/■*  K) 


2690 


530 


27800 


9280 
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